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ABSTRACT 


Infrared  polarization  techniques  to  improve  the  target-background  contrast  are  the 
subject  of  much  study  lately.  Polarized  infrared  images  of  the  research  vessel,  R/V 
POINT  SUR  were  taken  extensively  during  the  Electro-Optical  Propagation  Assessment 
in  Coastal  Environment  (EOPACE)  operational  period  in  Mar- Apr  1996  at  Point  Loma, 
San  Diego.  Contrast  improvement  analysis  of  these  images  requires  an  infrared  model  of 
the  research  vessel. 

This  thesis  models  the  apparent  infrared  signature  of  the  RfV  POINT  SUR  using 
ship  skin  temperature  records  from  the  EOPACE  measurements.  Using  a  simple 
geometric  model  and  adopting  a  pixel-by-pixel  approach,  the  emitted  radiation  for  each 
pixel  surface  is  calculated  from  the  Planck  radiation  law.  Radiation  reflected  off  the  pixel 
surface  is  estimated  using  the  polarized  SeaRad  atmospheric  propagation  code.  The  total 
radiance  is  compensated  for  atmospheric  transmittance  with  the  Navy  Aerosol  Model  to 
arrive  at  the  apparent  radiance  at  the  sensor.  The  apparent  radiance  is  translated  into 
apparent  temperature  for  comparison  with  the  recorded  images  to  validate  the  accuracy  of 
the  model.  There  is  moderate  agreement  between  the  model  and  recorded  images  with 
45-50%  of  the  pixels  falling  within  1.75K  of  the  measured  apparent  temperatures. 
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I.  INTRODUCTION 


Every  object  is  a  sotirce  of  infrared  radiation  with  no  exception  even  to  the 
atmosphere  through  which  the  radiation  propagates.  Therefore,  the  design  of  infrared 
detection  and  tracking  systems  requires  the  ability  to  discriminate  the  object  of  interest 
from  a  highly  cluttered  environment.  Various  filtering  techniques  such  as  spectral  and 
spatial  filtering  are  typically  used  to  achieve  a  detectable  signal  difference  between  the 
target  and  the  background.  More  recently,  polarization  techniques  to  improve  the  target  to 
background  contrast  in  the  infrared  wavelengths  are  being  seriously  studied.  Preliminary 
studies  done  at  the  Naval  Postgraduate  School  using  the  Marine  Aerosol  Properties  and 
Thermal  Imager  Performance  (MAPTIP)  trial  data  showed  that  the  ship  and  the  sky 
backgroimd  exhibit  no  significant  polarization  features  while  the  sea  exhibits  a 
considerable  degree  of  vertical  polarization  at  near  grazing  angles.  [Ref.l  and  Ref  2]  The 
studies  also  showed  that  horizontal  polarization  filtering  gives  a  10  to  20%  improvement 
for  ship/sea  contrast  and  up  to  15%  improvement  for  sea/sky  contrast  at  the  horizon.  The 
analysis,  however,  was  based  on  apparent  radiance,  that  is,  radiance  measured  at  the 
sensor.  It  does  not  show  how  the  relative  contribution  to  the  degree  of  polarization  from 
the  various  sources  and  atmosphere  add  up  to  give  the  observed  polarization  features.  A 
more  conclusive  study  would  need  to  accoimt  for  the  polarization  features  based  on 
actual  radiance  derived  from  radiation  theory.  Such  a  study  requires  good  environmental 
data  and  a  reliable  infrared  model  of  the  ship. 

Electro-optical  Propagation  Assessment  in  Coastal  Environment  (EOPACE) 
organized  by  NCCOSC-NRaD  (now,  SpaWarSysCen.San  Diego)  is  a  5-year  multi¬ 
national  effort  to  improve  the  performance  assessment  for  electro-optical  systems 
operating  in  a  coastal  environment.  This  is  because  the  coastal  environment  may  differ 
significantly  from  open  ocean  conditions,  and  existing  propagation  codes  such  as 
LOWTRAN/MODTRAN  may  not  adequately  account  for  the  variations.  While  the  main 
effort  of  the  EOPACE  measurement  series  is  to  fully  characterize  the  coastal 
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environment,  it  also  provided  the  opportunity  to  study  the  performance  of  polarization 
techniques.  During  one  of  the  EOPACE  operational  periods  in  Mar-April  1996  at  Point 
Loma,  San  Diego,  extensive  polarized  infrared  images  of  the  research  vessel,  RA^  POINT 
SUR  were  taken.  With  the  comprehensive  environmental  database  gathered  during  that 
period,  it  is  possible  to  perform  actual  radiance  calculations  for  the  ship  and  the  marine 
background. 

This  thesis  models  the  apparent  infrared  signature  of  the  research  vessel  RW 
POINT  SUR  using  ship  skin  temperature  records  from  the  EOPACE  measurements.  The 
research  vessel  is  represented  by  a  geometric  model  consisting  of  23  vertical  planes.  Each 
plane  is  associated  with  a  temperature  reading  from  one  of  the  14  thermistors  located 
around  the  lUV  POINT  SUR  during  the  imaging  period.  To  facilitate  subsequent 
comparison  with  the  recorded  images,  each  plane  is  broken  into  pixel  size  surfaces.  For 
each  pixel  size  surface,  the  emitted  radiation  from  the  surface  is  calculated  from  the 
Planck  radiation  law.  The  radiation  reflected  off  the  pixel  surface  into  the  sensor  is 
estimated  from  the  environmental  data  using  SeaRad  [Ref  3].  SeaRad  is  a  background 
add-on  to  the  LOWTRAN  7  /MODTRAN  atmospheric  codes  that  calculates  the  sea 
radiance  based  on  a  Cox  and  Munk  [Ref. 4]  wind  driven  capillary  wave  slope  distribution. 
In  the  polarized  SeaRad  version,  the  total  sea  radiance  is  broken  down  into  its  vertically 
and  horizontally  polarized  components.  This  feature  allows  us  to  derive  the  apparent 
radiance  separately  for  the  horizontal  and  vertical  polarizations.  For  atmospheric 
compensation,  the  Navy  Aerosol  Model  (NAM)  developed  specially  for  marine  boxmdary 
layer  propagation  is  used.  Finally,  the  model  is  reconstructed  into  a  2-dimensional  array 
with  the  elements  giving  the  apparent  radiance  values  for  each  pixel.  The  array  is  then 
compared  with  the  ship  pixels  of  the  image  arrays  for  each  polarization  to  validate  its 
accuracy.  A  reliable  mfrared  model  of  the  research  vessel  will  offer  many  more  options 
for  future  analysis,  such  as  contrast  improvement  analysis,  range  improvement  analysis 
and  sensitivity  analysis. 
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II.  MODELING  APPROACH 


A.  SURFACE  RADIOSITY  EQUATION 

Modeling  the  infrared  signature  of  a  target  or  background  is  a  complex  process 
because  every  object  is  a  source  of  both  emission  and  reflection  radiation.  Hence,  all 
infrared  models  are  concerned  with  attempting  to  solve  the  surface  radiosity  equation 
[Ref  5]  : 


where 


Ls  =  spectral  radiance  of  the  siirface  S, 

6r  =  zenith  angle  of  the  reflected  radiation  towards  the  sensor, 

(j)r  =  azimuth  angle  of  the  reflected  radiation  towards  the  sensor, 

0i  =  zenith  angle  of  the  incident  radiation, 

(j)!  =  azimuth  angle  of  incident  radiation, 

X  =  wavelength  of  the  radiation, 

Ss  ==  emissivity  of  the  surface, 

B  =  Planck 's  function, 

Ts  =  surface  temperature, 

fi  =  Bi-directional  Reflectance  Distribution  Function  (BRDF)  of  the  surface, 
Li  =  incident  spectral  radiance,  and 
Qi  =  incident  projected  solid  angle 


This  equation  essentially  states  that  the  spectral  radiance  of  a  surface  element,  S, 
is  the  sum  of  two  components;  the  emitted  radiation  and  the  reflected  radiation.  The 
emitted  radiation  is  the  product  of  the  Planck  function  at  the  temperature  of  the  surface 
and  the  directional  spectral  emissivity  of  the  surface.  The  reflected  radiation  is  the 
integral  over  all  the  incident  directions  of  the  product  of  the  incident  spectral  radiance  for 
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an  elemental  projected  solid  angle  of  the  source  and  the  surface’s  BRDF. 


dQj 


Sensor 


Figure  2.1.  Diagram  of  geometric  relationships  [from  Ref  5] 

The  surface  radiosity  equation  may  seem  to  be  a  relatively  simple  equation  to 
solve.  However,  this  is  not  the  case.  For  instance,  there  are  inter-body  reflections.  Inter¬ 
body  reflections  involve  radiation  that  is  reflected  to  and  fro  between  two  surfaces. 

Hence,  part  of  the  incident  radiance  onto  a  surface  may  be  radiance  from  the  sxirface  itself 
reflected  off  a  second  surface  back  to  itself  and  that  portion  of  the  incident  radiance  in 
imknown.  Furthermore,  performing  the  integral  for  reflected  radiance  is  not  easy  and 
usually  requires  numerical  methods.  Given  the  limited  time  and  resources  available  for 
this  thesis,  a  more  practical  approach  is  to  simplify  the  equation  to  a  more  useable  form 
with  a  few  assumptions. 

Firstly,  we  see  that  inter-body  reflections  are  at  least  twice  reflected.  In  most 
cases,  the  reflectivity  is  low  and  hence,  inter-body  reflections  can  be  reasonably 
neglected.  Secondly,  by  assuming  that  the  surface  is  perfectly  diffuse,  we  can  use  the 
Lambertian  approximation  to  eliminate  the  angular  dependence  of  the  emissivity  and 
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reflectivity.  Thirdly,  assuming  that  the  object  behaves  like  a  graybody,  and  making  use  of 
Kirchoff  s  law,  we  eliminate  spectral  dependence  of  the  emissivity  and  reflectivity. 

Hence,  we  obtain  a  simpler  equation : 


L,(X)  =  s.B,(T,-A)+^-e.)\L,(e,4,-,X) 


We  are  still  left  with  the  task  of  performing  the  integral  for  reflected  radiance, 
which  is  an  angular  average  of  the  incident  environmental  radiance  .  For  each  surface 
element  orientation,  this  integral  average  will  be  different  and  will  have  to  be  taken  into 
account  for  a  collection  of  elemental  siufaces  that  may  represent  an  object.  A  further 
simplification  is  needed  and  this  will  be  discussed  in  the  later  section  of  this  chapter. 


1.  Geometric  Representation 


To  apply  the  surface  radiosity  equation  to  the  problem ,  three  basic  issues  would 
have  to  be  settled  first.  The  first  issue  is  the  geometric  representation  of  the  object  of 
interest.  The  object  of  interest  is  the  135-feet  research  vessel,  RW  POESIT  SUR,  owned 
by  the  National  Science  Foimdation  and  operated  by  the  Central  California 
Oceanographic  Cooperative  at  Moss  Landing  Marine  Laboratories,  Monterey.  During  the 
recording  period,  it  was  used  as  the  cooperative  target  for  infrared  imaging.  The  general 
data  of  the  ship  is  given  in  Appendix  A. 


The  simplest  and  most  convenient  way  to  represent  any  object  is  to  use  a 
collection  of  flat  facets  wherever  possible.  In  fact,  it  is  true  to  say  that  this  is  done  in 
many  infrared  models.  As  observed  from  the  infrared  images  of  the  ship  at  ranges  of 
about  one  kilometer  and  more,  intricate  details  of  the  ship  are  not  readily  apparent. 
Furthermore,  given  the  relatively  small  look  down  angle  at  the  sensor,  all  horizontal 
planes  have  a  small  viewable  area,  and  do  not  contribute  significantly  to  the  emitted 
radiance.  However,  it  should  be  noted  that  at  grazing  angles,  radiance  reflected  off 
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horizontal  planes  may  be  significant.  Nevertheless,  we  will  neglect  all  horizontal  planes 
and  construct  a  simple  scaled  model  of  the  ship  using  a  collection  of  vertical  planes.  The 
geometric  model  is  shown  in  Figure  2.2a  and  2.2b  below.  The  mast  structure  of  the  ship 
is  not  included  as  it  is  barely  visible  in  most  of  the  thermal  images  in  this  study.  As 
mentioned  above,  the  details  at  the  bridge  and  01  deck  do  not  show  up  in  the  thermal 
images  and  they  are  modeled  as  simple  cuboids  here.  The  bow  section  is  approximated  as 
2  regular  vertical  planes  (planes  22  and  23).  This  makes  the  geometry  analysis  for  the 
bow  section  much  simpler.  There  are  a  total  of  23  vertical  planes,  numbered  as  shown  in  . 
the  figures.  The  dimensions  of  each  plane  are  given  in  Appendix  A  as  well. 


Figure  2.2a.  Port  view  of  the  geometric  model,  and  the  plane-thermistor  matching  of 
RA^  POINT  SUR.  The  index  number  such  as  "Til"  refers  to  the  thermistor  number  as 
listed  in  Appendix  B. 


6 


Figure  2.2b.  Starboard  view  of  the  geometric  model  and  the  plane-thermistor  matching  of 
RA^  POINT  SUR 

2.  Surface  Temperature 

The  second  issue  is  to  decide  how  these  planes  obtain  their  temperatures.  This 
will  allow  us  to  calculate  the  emission  term  in  the  surface  radiosity  equation.  One  would 
have  to  consider  the  sources  of  heat  generation  on  the  ship,  such  as  engines,  running 
machineries,  electrical  switchboards,  etc.,  the  various  heat  transfer  processes  and  also 
interaction  with  the  environment  to  arrive  at  a  "first  principles"  temperature  distribution 
for  the  ship  surfaces.  This  is  by  no  means  a  simple  task  and  is  done  in  more  advanced 
models  such  as  Target  Contrast  Model  2  (TCM2)  [Ref.6]  adapted  for  use  in  the  EOTDA 
FLIR  Tactical  Decision  Aid  Program.  For  our  purposes,  we  require  a  simpler  model. 
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For  the  period  of  the  EOF  ACE  imaging  measurements,  14  thermistors  were 
attached  at  selected  locations  around  the  RA^  POINT  SUR.  The  exact  locations  of  the 
thermistors  are  listed  in  Appendix  B.  The  thermistor  readings  were  recorded  at  20-second 
intervals  throughout  the  whole  imaging  period.  Hence,  we  can  associate  each  plane  of  our 
geometric  model  with  a  temperature  from  one  of  the  14  thermistor  readings. 

In  associating  the  planes  to  the  temperature  reading  of  a  thermistor,  it  is  also 
helpful  to  compare  the  apparent  temperature  distribution  from  the  thermal  images  to  the 
thermistor  readings.  Figure  2.3  below  shows  a  sample  thermal  image.  It  can  be  seen  from 
the  image  that  there  is  a  hot  section  that  extends  from  the  stack  down  to  the  hull  on  the 
port  side  of  the  ship.  As  a  result,  an  additional  plane  was  inserted  (plane  2)  to  better 
reflect  the  temperature  distribution  of  the  hull  on  the  port  side.  (See  figure  2.2a)  The 
matching  of  each  plane  to  the  temperature  reading  of  a  thermistor  is  indicated  beside  the 
plane  number  in  Figure  2.2a  and  2.2b  above. 


Figure  2.3  A  typical  portside  view  thermal  image  from  file  "Based  .pte".  A  hot 
section  extending  from  the  stack  down  to  the  hull  of  the  ship  is  visible 


3.  Radiation  from  the  Environment 

There  are  two  contributors  to  the  reflected  radiation  term  in  the  surface  radiosity 
equation.  Since  we  have  ignored  inter-body  reflections,  the  environment  is  the  sole 
source  of  reflected  radiation. 
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Modeling  of  the  environment  is  yet  in  itself  another  complex  process.  SeaRad 
with  the  Navy  Aerosol  Model  is  particularly  suited  for  this  case  since  it  calculates  sea 
radiance,  reflected  sky  radiance,  solar  radiance,  atmospheric  path  radiance,  scattering  and 
transmittance  over  the  marine  boundary  layer  atmosphere.  However,  like  LOWTRAN  7 
and  MODTRAN,  each  input  case  is  restricted  to  a  point  to  point  calculation.  What  this 
means  is  that  the  solid  angle  subtended  by  the  source  area  has  to  be  reasonably  small  for 
that  calculation  to  be  valid.  Therefore,  the  correct  approach  is  to  consider  the  reflection 
off  an  elemental  surface  for  each  plane.  To  facilitate  the  subsequent  comparison  with  the 
recorded  images,  the  size  of  each  elemental  surface  is  made  to  coincide  with  the  pixel 
size  of  the  image. 

Going  back  to  our  simplified  equation  2.2,  we  are  still  left  with  the  task  of 
performing  the  integral  for  reflected  radiance  over  the  hemisphere  of  the  elemental 
surface  exposed  to  the  environment.  For  each  elemental  surface,  this  integral  would 
require  approximately  50,000  input  cases  to  SeaRad.  Given  the  dimensions  of  RW 
POINT  SUR  and  the  imaging  ranges,  it  is  estimated  that  each  frame  will  have  an  average- 
of  250  ship  pixels.  Thus,  a  huge  amount  of  computational  resources  is  needed  for  proper 
modeling.  A  ftirther  simplification  is  required  here. 

Consider  a  pixel  surface,  which  is  at  a  height,  h,  above  the  sea  level.  The  range  to 
the  horizon  for  the  given  height  h  of  the  pixel  is  V(Rh)  where  R  is  the  radius  of  the  earth. 
The  angle  of  elevation  of  the  horizon  at  the  pixel  is  then  given  by  tan'*  VCR/h).  (See  figure 
2.4  below)  Since  h  is  very  small  compared  to  R,  the  value  of  V(R/h)  is  very  large.  This 
gives  a  to  be  nearly  90  degrees.  Therefore,  each  pixel  will  see  the  sky  and  sea  as 
approximately  two  equal  half  hemispheres  each  subtending  a  solid  angle  of  n  steradians. 
If  we  assume  that  the  incident  sea  and  sky  radiance  are  constant  for  their  respective  half 
hemispheres,  then  equation  2.2  becomes  : 


V 


(2.3) 
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where  Lsea  is  the  average  sea  radiance  over  the  lower  half  hemisphere  and 
Lsky  is  the  average  sky  radiance  over  the  upper  half  hemisphere. 


Figure  2.4  Geometry  relation  of  a  pixel  with  the  sky  and  sea.  Each  pixel  will  see  the 
horizon  at  nearly  90  degrees  and  hence,  the  sea  and  sky  as  two  rougly  equal  half 
hemispheres. 


For  a  cloudless  sky,  the  sky  radiance  may  be  considered  to  be  the  sum  of  the  sky 
background  emission  and  the  pafli  radiance.  It  is  reasonable  to  assiune  that  the  sky 
background  is  isotropic.  Therefore,  sky  background  emission  is  independent  of  look 
angle.  However,  the  path  length  and  hence,  the  atmospheric  path  radiance  would  increase 
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with  increasing  zenith  angle  towards  the  horizon.  Hence,  the  most  appropriate  path  to  be 
used  for  SeaRad  calculation  so  as  to  obtain  the  average  sky  radiance  should  be  the  path 
that  gives  the  average  path  radiance.  This  was  found  to  be  the  slant  path  to  space  at  a 
zenith  angle  of  65.3  degrees. 

For  the  average  sea  radiance,  we  will  assume  the  sea  to  be  a  flat  calm  surface.  The 
sea  radiance  is  once  again  given  by  the  surface  radiosity  equation  2. 1 .  We  can  write  the 
simplified  equation  as  : 


Lsea  Bsea "  r(oc)(B  sea  "  Lsky) 


(2.4) 


where  r  is  the  reflectance  given  by  the  Fresnel  equations , 
a  is  the  angle  of  incidence  or  elevation  and 
we  have  assumed  here  that  Lsky  is  the  average  sky  radiance. 

We  can  see  that  Lsea  is  also  a  function  of  the  zenidi  angle.  Since  Bsea  and  Lsky  are  both 
constant,  the  average  Lsea  is  then  given  by  the  path  for  average  reflectance.  The  path  that 
would  best  represent  the  average  sea  radiance  was  found  to  be  the  slant  path  to  the  sea  at 
a  zenith  angle  of  1 18.1  degree.  The  details  of  both  approximations  are  given  in  Appendix 
C. 
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in.  IMPLEMENTATION 


The  modeling  and  analysis  programs  for  this  thesis  were  written  in  Interactive 
Data  Language  (IDL).  However,  SeaRad  is  a  FORTRAN  code  written  to  run  in  an  MS- 
DOS  environment.  Therefore,  the  implementation  stage  has  to  be  divided  into  3  parts, 
switching  from  IDL  to  MS-DOS  and  then  back  to  the  IDL  environment. 

In  the  first  part,  the  programs  perform  3  main  functions.  These  are  : 

-  to  extract  the  relevant  meteorological  data  required  for  SeaRad, 

-  to  calculate  the  path  geometry  between  each  pixel  and  the  sensor,  and 

-  to  write  the  meteorological  and  geometry  data  into  the  required  input  format  for 
SeaRad. 

For  the  second  part,  the  main  task  is  to  run  the  input  files  through  the  SeaRad 
code.  After  running  SeaRad,  the  tasks  of  the  third  part  are  : 

-  to  extract  the  required  data  from  the  SeaRad  output  files 

-  to  calculate  the  actual  radiance  received  at  the  sensor 

-  to  build  a  2  dimensional  array  of  the  model  that  v^dll  facilitate  subsequent 
comparison  and  analysis. 

A  list  of  the  programs  written  for  each  part  is  given  in  Appendix  D.  The  details  of 
operations  and  information  for  each  part  are  discussed  in  the  following  sections. 

A.  PART  1  -  EXTRACT  METEOROLOGICAL  DATA,  CALCULATE  PATH 
GEOMETRY  AND  WRITE  SEARAD  INPUT  FILES 

There  are  4  main  procedures  written  for  this  part  and  they  run  in  the  sequence  ; 
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START.PRO,  PLANEIN.PRO,  SEAIN.PRO,  AVESEAIN.PRO.  Each  main  procedure 
performs  a  specific  function.  START.PRO  extracts  the  relevant  meteorological  dgta  It 
also  provides  the  basic  geometry  parameters  to  the  procedure  PLANEIN.PRO,  which 
calculates  the  pixel  coordinates  and  pixel-sensor  path  geometry.  After  all  the  required 
data  are  extracted  and  calculated,  the  procedures  SEAIN.PRO  and  AVESEAIN.PRO  are 
used  to  write  the  data  into  FORTRAN  input  format  required  to  nm  the  SeaRad  code.  The 
details  given  in  the  subsections  below  will  more  or  less  follow  this  sequence  of  operation. 

1.  Data  File  Format 

The  original  infrared  images  of  the  RA^  POINT  SUR  were  recorded  using  the 
PTRWIN,  version  3.16  data  acquisition  software  from  CEDIP.  They  are  identified  with 
the  file  extension  .PTW.  Although  PTRWIN  has  the  capability  to  analyze  the  images,  it  is 
inadequate  for  analysis  with  large  numbers  of  images.  Therefore,  they  were  converted 
into  a  more  efficient  format  using  an  IDL  procedure  written  by  Pontes  [Ref  7]  called 
BASEGEN.PRO.  The  converted  files,  identified  by  the  filename  "BASE*  .PTE",  have  a 
basic  structure  that  contains  an  information  vector,  a  time  vector  that  gives  the  time  of 
each  frame  and  is  followed  by  the  images  in  the  corresponding  order.  The  detailed 
structure  and  type  of  information  carried  by  each  "BASE*  .PTE"  file  is  given  in  Appendix 
E. 


Each  "BASE*.PTE"  carries  a  total  of  approximately  75  picture  frames,  that  is, 
roughly  25  frames  for  each  case  of  horizontally  polarized,  vertically  polarized  and 
unpolarized  images.  Given  that  the  scaimer  produces  6.25  picture  frames  per  second  and 
that  each  image  set  of  different  polarization  is  taken  sequentially,  there  is  a  time  lag  of 
about  8  seconds  between  the  first  frame  and  the  last  frame  of  the  file.  One  of  the 
objectives  in  the  previous  analysis  performed  by  Pontes  [Ref  7]  on  the  EOPACE  data 
was  to  determine  whether  there  is  a  ship  radiance  variation  with  time.  It  was  found  in  that 
analysis  that  the  variation  of  radiance  within  the  time  scale  of  8  seconds  is  not  significant. 
Therefore,  only  the  first  frame  of  each  polarization  case  will  be  considered  in  this  thesis. 
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2.  Geometry  of  Measurements 


The  measurements  were  made  from  the  NCCOSC  Building  15  (32°  39'  55"  N, 

117°  14'  32"  W,  elevation  above  sea  level  -57  feet)  at  the  tip  of  Point  Loma.  The  research 
vessel  was  imaged  at  different  aspect  angles  at  ranges  from  0.5  km  up  to  5  km.  The 
bearing  of  the  ship  was  also  made  to  vary  from  about  100  to  260  degrees  from  the 
observation  point.  A  chart  showing  the  location  of  the  measurement  is  in  Appendix  F. 

3.  Meteorological  Data 

Meteorological  data  were  collected  at  2  separate  meteorological  towers  (METOC 

1  and  METOC  2)  onboard  the  R/V  POINT  SUR.  The  data  from  METOC  1  and  METOC 

2  were  recorded  in  files  identified  by  "EOP*.TXT"  and  "EOP30SEC.*"  respectively, 
where  *  indicates  the  day  of  the  year.  Each  file  was  started  at  0:00  Greenwich  Mean 
Time  (GMT)  of  each  GMT  day  and  the  data  were  recorded  in  30-second  intervals.  For 
each  30-second  interval,  23  parameters  were  recorded.  These  are  the  GMT  time  of  the 
reading,  relative  wind  speed  and  direction,  ship  speed  and  direction,  true  wind  speed  and 
direction,  air  temperature,  relative  humidity,  pressure,  sea  surface  temperature  and  the 
GPS  position  of  the  ship.  The  detailed  format  of  the  METOC  files  is  included  as 
Appendix  G. 

To  extract  the  meteorological  data,  the  time  vector  in  the  .PTE  file  is  first  read. 
This  is  then  converted  to  GMT  to  open  the  correct  files  and  to  read  the  nearest  line  of 
data  recorded  for  that  time.  Since  not  all  the  data  are  needed,  the  relevant  data  are 
extracted,  averaged  and  used  to  calculate  some  additional  required  information.  The 
program  READMET.PRO  does  all  these,  and  returns  a  vector  containing  the  range  of 
ship,  sea  temperature,  air  mass  parameter,  current  wind  speed,  24  hour  wind  speed,  wind 
azimuth  relative  to  observer  azimuth,  GMT  day  of  the  year,  GMT  and  observer  azimuth. 
This  information  for  the  analyzed  files  is  documented  in  Appendix  H. 
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The  US  Navy  Aerosol  Model  (NAM)  was  developed  to  compute  the  aerosol 
extinction  coefficients  for  visual  and  infrared  wavelengths  in  the  marine  boundary  layer. 
Since  we  are  concerned  with  the  propagation  of  infrared  radiation  from  a  ship  in  a  coastal 
environment ,  it  is  particularly  suited  for  application  in  this  thesis.  However,  in  order  to 
use  the  Navy  aerosol  model,  part  of  the  information  required  is  the  "air  mass  parameter 
(AMP)",  which  is  included  in  the  MODTRAN  propagation  code  as  "ICSTL".  This 
parameter  is  not  recorded  in  the  METOC  files  and  cannot  be  derived  from  any  of  the 
parameters  collected  there. 


The  traditional  way  to  determine  the  air  mass  parameter  is  to  measure  the  gas 
concentration  of  the  radioactive  Radon-222  in  the  atmosphere.  During  the  Intensive 
Operational  Period  (lOP)  in  Mar- April  1996,  the  Radon-222  gas  concentration  was 
measured  at  San  Nicolas  Island,  Point  Mugu  and  San  Diego.  This  information  is  available 
to  the  EOP ACE  participants  at  the  EOPACE  website  at  http://sunspot.spawar.navy.mil:  80 
/ 543/ eopace/ sd_iop_ma96/nrad883/bdg323/radon/radonsd.txt" .  With  these 
measurements,  the  air  mass  parameter  is  then  simply  given  by  the  equation  [Ref  9]  : 


P  =  INT\ 


4 


+  1 


(3.1) 


where  INT  truncates  to  the  integer  value  and  Rn  is  the  atmospheric  Radon-222  gas 
concentration  in  picocurries  per  cm^.  It  should  be  noted  here  that  this  equation  is 
empirical  and  would  have  to  be  modified  for  application  in  other  geographical  regions.  A 
brief  discussion  on  the  NAM  and  the  air  mass  parameter  is  in  Appendix  I. 

4.  Pixel  Size 

The  infrared  scanner  used  to  record  the  ship  images  is  the  dual  band  serial 
scanning  AGA-780  Thermovision  Thermal  Imaging  System.  The  short-wave  (SW) 
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scanner  uses  an  indium  antimonide  (InSb)  detector  and  covers  the  3-5.6  )j.m  spectral 
band.  The  long- wave  (LW)  scaimer  covers  the  8-14  pm  spectral  band  and  has  a  mercury- 
cadmixrm-telluride  (HgCdTe)  detector. 


Electro-magnetic  radiation  from  an  object  is  focused  onto  the  detectors  by  two 
separate  systems  of  infrared  lenses  and  rotating  prisrhs.  In  each  scanner,  an  infrared 
image  frame  is  produced  by  2  prisms,  rotating  in  such  a  way  that  they  produce  four  fields 
to  make  one  interlaced  frame.  This  is  shown  in  Figure.3.1  below.  Each  field  has  100 
horizontal  scanning  lines  and  only  70  of  these  are  used  as  active  imaging  lines.  This 
gives  a  vertical  scanning  efficiency  (rise)  of  70%. 


Figure  3.1  Interlace  fields  of  the  AGA-780  [Ref.9] 

There  were  3  different  lenses  that  were  available  for  use  with  the  scanner.  These 
are  the  3.5°,  7°  and  20°  field  of  view  (FOV)  lenses.  Table  3.1  lists  the  specifications  for 
these  lenses. 
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Table  3.1.  Lens  specification  from  AGA  780  Manual  [Ref.9] 


Lens 

Min  focus 

Focal  distance 

Geometric  resolution 

3.5° 

3.3m 

191mm 

0.5mrad 

7° 

1.0m 

99mm 

l.lmrad 

20° 

0.5m 

33mm 

3.4mrad 

For  the  whole  recording  period  during  EOF  ACE  April  96,  only  the  3.5°  and  7° 
lenses  were  used.  The  horizontal  pixel  resolution  is  therefore  O.Smrad  and  1.1  mrad 
respectively.  To  obtain  the  effective  vertical  pixel  resolution,  we  multiply  the  geometric 
resolution  by  the  scanning  efficiency  (t|sc)  .  This  gives  0.35mrad  and  0.77mrad  for  the 
3.5°  and  7°  FOV  respectively. 

5.  Path  Geometry 

The  idea  behind  the  path  geometry  calculation  was  to  first  determine  the 
coordinates  of  the  center  of  each  pixel  (which  are  called  nodes)  with  respect  to  an  origin 
vertically  below  the  sensor.  Thence,  the  sensor  to  pixel  path  is  a  vector  given  by  the 
difference  of  the  pixel  coordinates  and  sensor  coordinates. 

To  obtain  the  3  dimensional  coordinates  of  each  node,  we  need  to  consider  each 
plane  of  the  geometric  model  of  the  ship  (Figure  2.2,  Chapter  II)  one  at  a  time.  The 
program  PLANE.PRO  is  first  used  to  select  the  visible  planes  for  computation.  This 
significantly  reduces  the  computation  load.  Each  selected  plane  is  then  broken  down  into 
nodes  depending  on  the  size  of  the  elemental  surface,  which  in  turn,  is  a  fimction  of  the 
range  and  pixel  resolution.  The  position  of  each  node  relative  to  the  center  of  the  ship  is 
first  calculated.  A  rotation  is  then  applied  to  the  node  positions  to  orientate  the  plane  to 
the  correct  aspect  angle.  Lastly,  the  origin  is  shifted  from  the  center  of  the  ship  to  an 
origin  directly  below  the  sensor. 
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One  of  the  main  problems  at  this  stage  is  the  fractional  edges  of  the  planes.  As 
each  plane  may  not  contain  whole  number  of  pixels,  the  edges  have  to  be  carefully 
handled  so  that  when  the  neighboring  plane  is  put  together,  it  does  not  give  an  incorrect 
radiance  value  later  on.  To  do  this,  the  pixels  at  the  edge  are  always  rounded  up.  Then  a 
weighting  factor  based  on  remaining  area  at  the  edge  to  the  pixel  area  is  calculated.  For 
those  pixels  at  the  edge  of  a  plane,  the  weighting  would  be  less  than  one  while  the 
weighting  for  the  pixels  at  the  edge  of  the  adjacent  planes  will  be  one  minus  the 
weighting  of  the  edge  pixels  of  the  first  plane.  These  pixels  will  share  the  same  node 
positions.  Therefore,  the  ouQ)Ut  of  the  routine  PLANEIN.PRO  contains  a  weighting 
matrix  in  addition  to  the  3  dimensional  coordinates  of  the  position  of  the  nodes. 

The  reflection  path  was  discussed  towards  the  end  of  Chapter  11.  The  paths 
representing  the  average  sea  and  sky  are  just  the  slant  paths  with  zenith  angles  of  1 1 8. 1 
and  65.3  degrees  at  the  pixels  respectively. 

6.  SeaRad  [Ref.  3] 

SeaRad  was  developed  as  an  ocean  background  add-on  to  the  MODTRAN  2  and 
LOWTRAN  7  atmospheric  codes.  The  code  used  to  calculate  atmospheric  transmittance, 
solar,  sky  and  path  radiance  is  the  same  code  as  MODTRAN/LOWTRAN.  The  sea 
radiance  add-on  is  activated  by  a  "seaswitch".  When  this  switch  is  set  to  true,  the  code 
will  take  the  boundary  temperature  as  the  sea  surface  temperature  whenever  the  path 
intersects  the  sea  surface.  It  then  applies  the  Fresnel  coefficients  to  the  reflection  of 
infrared  radiation  by  a  seawater  capillary  wave  facet.  The  result  of  a  single  tilted  facet  is 
then  integrated  over  all  possible  tilts  which  follows  a  Cox  and  Munk  probability 
distribution  function.  Hence,  SeaRad  calculation  gives  a  pixel  average  sea  radiance, 
which  suited  our  purpose  here. 

In  addition,  it  is  worth  noting  that  SeaRad  does  not  take  into  account  the 
correction  for  shadowing  effects  proposed  by  Saunders  [Ref.  10].  In  essence,  Saunders 
reasoned  that  radiation  from  certain  portions  of  a  wave  is  blocked  by  the  wave  in  front. 
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Hence,  these  portions  do  not  contribute  to  the  sea  radiance.  The  amount  of  area  being 
blocked  is  dependent  on  the  viewing  angle.  The  more  acute  the  viewing  angle,  the  more 
area  of  the  wave  is  blocked.  Given  the  viewing  geometry  of  our  problem,  this  effect 
maybe  appreciable. 

SeaRad  input  requires  the  submission  of  a  series  of  at  least  5  cards.  Each  card  is  a 
line  of  formatted  input.  As  each  pixel  requires  one  set  of  input  cards  and  there  is  a  large 
number  of  pixels,  the  routines  SEAIN.PRO  and  AVESEAIN.PRO  are  used  to  write  all 
the  extracted  meteorological  data  and  calculated  geometry  data  into  the  required  format 
for  each  pixel.  A  sample  input  file  is  shown  at  Appendix  J.  An  explanation  of  the  input 
cards  used  for  the  input  files  is  also  included  in  Appendix  J. 

B.  PART  2  -  RUN  SEARAD 

SeaRad  is  run  in  an  MS-DOS  environment.  The  DOS  commands  to  run  a  single 
input  file  are  : 

copy  "input  filename"  tape  5 

SeaRad 

The  output  of  the  calculations  is  written  to  an  output  file  called  OUT.  As  can  be 
seen,  SeaRad  is  not  designed  to  run  a  large  number  of  input  files  as  is  needed  here.  Since 
each  input  file  differs  only  in  the  geometry  data,  it  would  be  possible  to  modify  the 
SeaRad  codes  to  accept  multiple  geometry  inputs.  A  simpler  alternative,  though  less 
efficient,  is  to  run  SeaRad  in  batches  using  MS-DOS  batch  file  commands.  A  batch  file 
called  SEAl  .BAT  was  written  to  run  input  files  in  batches.  There  are  3  types  of  input 
files  in  the  "CASEl/INPUTS/"  directory  :  the  first  type  is  for  the  pixel-sensor  path 
geometry ,  the  second  type  for  the  pixel-sea  geometry  and  the  last  type  for  the  pixel-sky 
geometry.  These  files  are  identified  by  the  filename  "*p*c**.std",  "*a*b*0.std"  and 
"*a*b*l.std"  respectively.  The  output  of  each  file  has  exactly  the  same  name  as  the  input 

file  but  they  are  copied  into  the  "  CASE  1 /OUTPUTS/"  directory. 
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C.  PART  3  -  EXTRACT  SEARAD  OUTPUT,  CALCULATE  RECEIVED 
RADIANCE  AND  CONSTRUCT  2-D  ARRAY 


After  running  SeaRad,  the  procedure  CONTINUE.PRO  reads  the  SeaRad  output 
files  to  extract  the  required  information.  It  then  calculates  the  radiance  received  at  the 
sensor  and  reconstructs  the  model  into  a  2  dimensional  array.  The  second  section  of 
CONTINUE.PRO  compares  the  resultant  array  with  the  recorded  image  arrays.  The 
second  section  will  be  discussed  in  the  next  chapter. 

1.  SeaRad  Output 

The  format  of  the  SeaRad  output  files  depend  on  the  type  of  calculation  that  is 
done.  For  the  three  input  cases,  two  types  of  output  files  were  written.  A  sample  of  each 
type  of  output  file  is  shown  in  Appendix  K.  The  routines  SEAOUT.PRO  and 
AVSEAOUT.PRO  are  used  to  extract  the  required  information  from  the  ou^ut  files.  The 
information  returned  by  these  programs  is  : 

-  pixel  to  sensor  path  transmittance  ( Xpath ), 

-  pixel  to  sensor  path  radiance  ( Lpath), 

-  total  radiance  at  the  pixel  for  the  average  sky  radiance  path  (Lsky) 

-  total  radiance  at  the  pixel  for  the  average  sea  radiance  path  in  horizontally 

polarized  (LHsea)-,  vertically  polarized  (Lvsea)-and  randomly  polarized  components 

(Lusea)- 


2.  Ship  Skin  Temperature 

As  previously  mentioned,  the  ship  skin  temperatures  of  RA^  POINT  SUR  were 
recorded  continuously  in  20-second  intervals  throughout  the  whole  measurement  period. 
These  readings  were  stored  in  a  file  called  "SKINDATA.PTE".  Each  line  has  16 
readings  corresponding  to  2  reference  temperature  readings  and  14  thermistor  readings. 
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The  detailed  format  for  the  skin  temperature  file  is  in  Appendix  L. 

3.  Ship  Skin  Radiance 

Given  the  ship  skin  temperature,  the  spectral  radiant  exitance  from  the  ship  skin 
surface  can  be  calculated  using  Planck  radiation  law  : 

<3.2) 

where  T  is  the  surface  temperature  in  K, 

X  is  the  wavelength  in  m 
h  is  Planck's  constant  (6.6256  x  lO'^'*  Js), 
c  is  the  speed  of  light  (3x10*  m/s), 
k  is  Boltzman's  constant  (1.38054  x  10'^*  J/K)  and 

Mbb  has  dimensions  in  watts  per  meter  square  per  wavelength  interval  in  m. 

The  radiant  exitance  for  a  graybody  over  %  steradians  of  solid  angle  is  then  given 
by  the  Stefan-Boltzman  law,  which  is  : 

Mcrayto<fyiT)  =  f  f  M ,,iA,T)dA  =  SOT^ 

where  s  is  the  emissivity  of  the  surface  and 

CT  is  the  Stefan-Boltzman  constant  (5.6697  x  10'^^  W/cm^K^). 

However,  we  are  only  interested  in  the  in-band  radiance  for  the  long-wave  band 
(8-12  pm).  The  in-band  radiance  can  be  represented  as  : 


^surface  =  “  I  T)dA 

7t 


(3.4) 
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There  are  two  methods  to  calculate  the  in-band  radiance.  One  method  is  to  use  the 
universal  blackbody  curve  or  pre-tabulated  integrals  over  it.  The  other  method  is  to 
perform  the  integration  numerically.  The  routine  "EMISSION.PRO"  performs  the 
numerical  integration  using  step  size  of  0.01  pm. 

Since  we  have  already  associated  each  plane  with  the  temperature  reading  of  a 
thermistor,  the  in-band  radiance  is  calculated  for  each  pixel  using  the  above  equation.  In 
the  calculations,  the  ship  surface  is  assumed  to  behave  like  a  graybody  with  an  emissivity 
previously  measured  to  be  0.95.  Since  no  further  information  on  the  optical  properties  of 
the  ship  paint  is  available,  another  assumption  was  made  that  the  horizontal  and  vertical 
polarized  components  of  the  in-barid  radiance  are  equal. 

4.  Radiance  Received  at  Sensor 

Finally,  the  radiance  received  at  the  sensor  is  the  sum  of  the  ship  skin  radiance 
compensated  for  atmospheric  transmittance,  the  path  radiance  and  the  reflected  sea  and 
sky  radiance.  This  is  best  summarized  by  the  following  formulae  for  the  two  polarized 
and  randomly  polarized  components  : 


L  =  — 

■*^Hsensor  ^ 


^ path^ surface  ^path 


path 


^  L  ^ 

L 

^Hsea  ^  ^ 

V  ^  J 


(3.5a) 


“*Vsensor 


path^ surface  ^path 


path 


^Vsea 


Z.  ^ 

si^ 


(3.5b) 


^Usensor  ~~  palh^ surface  ^path  "*■  ^  (3.5c) 

It  should  be  noted  here  that  there  would  be  a  transmission  loss  due  to  the  polarizing 
filters.  However,  it  is  not  necessary  to  correct  for  it  since  the  corrections  are  built  into  the 
scanner  calibration  process. 
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5.  Reconstructing  a  2-D  Array 


Remembering  that  the  procedure  was  done  one  plane  at  a  time,  the  remaining  task 
is  to  put  all  the  planes  together  to  give  a  2  dimensional  image  of  the  model.  This  process 
is  done  in  3  steps. 

The  first  step  is  to  multiply  the  radiance  matrix  of  each  plane  by  the  weighting 
matrix  so  that  the  pixels  at  the  edges  of  adjacent  planes  will  fit  together.  The  second  step 
is  to  correctly  position  each  plane  relative  to  a  reference  point  (the  origin  at  the  midpoint 
of  the  ship)  in  an  empty  72  x  30  matrix.  This  is  done  by  the  procedure  PLACE.PRO.  An 
additional  procedure  called  BLOCK.PRO  is  then  applied  to  blank  out  the  areas  of  overlap 
between  planes  for  that  viewing  angle.  Lastly,  the  radiance  values  of  the  planes  are 
simply  added  together  to  form  a  2D  array  of  ship  radiance  received  at  the  scanner.  A 
sample  image(magnified)  of  the  ship  model  derived  for  file  "BaseOl.pte"  is  shown  in 
Figure  3.2  below.  The  recorded  image  is  the  original  image  from  the  data  files.  The 
extracted  image  is  the  resultant  image  when  a  mask  having  the  same  shape  as  the  model 
but  containing  ones  and  zeros  is  applied  to  the  original  images. 


Recorded  image 


Model  Extracted  image 


Figure  3.2  Images  of  RA^  POINT  SUR  for  unpolarized  case  from  file  BaseOl.pte.  The 
model  image  shows  a  higher  temperature  than  the  recorded  temperature  in  the  extracted 
image.  Notice  that  some  of  the  background  pixels  are  inevitably  included  in  the  extracted 
image 
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IV.  COMPARISON  OF  RESULTS 


A.  APPROACH  FOR  COMPARISON 

1.  Pixel-by-pixel  Comparison 

Since  we  started  out  using  a  pixel  by  pixel  approach  to  derive  a  2-dimensional 
array  for  received  radiance  at  the  sensor,  each  element  in  the  2-D  array  will  directly 
correspond  to  a  pixel  in  the  recorded  image  frame.  Each  image  frame  has  a  sea  and  sky 
backgroimd.  Since  we  are  only  interested  in  comparing  the  ship  model  with  the  measured 
ship  radiance  the  background  would  have  to  be  eliminated.  To  extract  the  ship  pixels  in 
the  image  frame,  a  mask  is  first  created  from  the  70  by  30  model  array.  This  is  done  by 
setting  elements  of  ship  pixels  in  the  array  to  ones  and  the  backgroimd  to  zeros.  Next,  the 
coordinates  of  the  right  bottom  comer  of  the  ship  in  the  image  frame  are  read  from  user’s 
input  This  information  is  used  as  the  reference  point  for  aligning  the  mask  over  the 
image  frame  to  extract  the  ship  pixels  in  the  image.  Having  set  all  the  background 
elements  of  the  image  array  to  zero,  we  can  now  compare  the  two  arrays  directly  by 
simply  taking  their  difference.  However,  before  the  difference  can  be  taken,  the  elements 
in  the  arrays  must  have  the  same  units.  The  model  array  has  units  of  radiance  while  the 
image  array  has  units  of  digital  level.  A  unit  conversion  is  required. 

2.  Thermal  Level 

The  numerical  measure  of  the  received  and  detected  radiation  used  in  the  infrared 
scanner  is  called  the  thermal  level.  The  traditional  imit  for  the  thermal  level  is  the 
isotherm  unit,  which  is  linearly  related  to  the  received  photon  radiation.  Therefore,  the 
relationship  between  the  thermal  level  and  the  object  temperature  is  non-linear.  To 
translate  object  temperature  to  its  thermal  level,  a  calibration  function  is  used.  This 
calibration  function  is  obtained  by  plotting  the  thermal  level  against  the  known 
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temperature  of  a  blackbody.  It  has  a  characteristic  equation  of  a  form  similar  to  the 
Planck  equation : 


TL  =  -^ - OS  (4.1) 

-I 

where  A,  B,  C  are  constant  determined  by  curve  fitting, 

OS  is  the  offset  correction  for  detector  drift  and 
T  is  the  temperature  in  K. 

The  calibration  function  is  dependent  on  the  aperture  and  filter  used  and  a  set  of 
calibration  curves  is  usually  produced.  Prior  to  the  measurement  period,  the  scarmer  with 
polarizing  filters  was  calibrated.  The  calibration  curves  are  shown  in  Appendix  M.  The 
calibration  constants  for  each  polarization  filter  are  also  given  in  the  appendix.  Although 
all  the  calibration  constants  are  given,  it  should  be  noted  here  that  the  offset  constant  for 
detector  drift,  OS,  tends  to  vary  rapidly  with  changing  ambient  conditions.  Using  the 
calibration  functions,  the  thermal  level  of  the  ship  pixels  for  each  polarization  case  in  the 
images  can  be  easily  translated  to  the  equivalent  blackbody  temperature  of  the  ship 
surfaces. 


3.  Digital  Level 

It  was  found  that  the  images  recorded  by  the  PTRWIN  software  were  stored  in 
digital  levels  in  the  original  data  files.  Although  the  PTRWIN  software  have  the  function 
to  convert  digital  levels  to  isothermal  imits,  the  process  to  convert  all  the  "*.PTE"  files 
into  thermal  levels  would  require  re-storing  the  original  "*.PTW"  file  into  isothermal 
units  and  re-generating  the  "*.PTE"  files  using  BASEGEN.PRO.  An  alternative  way  is  to 
deduce  the  relationship  between  the  digital  level  and  thermal  level  from  the  PTRWIN 
software.  It  was  foimd  that  they  are  related  by  a  simple  linear  relation  which  depends  on 
the  Range  (R)  and  Level  (L)  settings  in  the  PTRWIN  software.  Hence,  to  convert  digital 
level  (DL)  to  thermal  level  (TL),  the  following  empirical  formula  is  used : 
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TL  =  A*DL  +  B 


(4.2) 


where  A  =  0.00122,  B-L-2.5  forR  — 5, 

A  =  0.00244,  B  =  L  -  5  for  R  =  10  and 

A  =  0.00488,  B  =  L-10  forR  =20. 

The  range  and  level  settings  are  not  available  in  the  information  vector  of  data 
(*.PTE)  files.  They  have  to  be  extracted  from  the  original  data  (*.PTW)  files.  Appendix 
N  documents  the  original  files  and  their  range/level  settings. 

3.  Calibration  Drifts 

The  time  lag  between  the  first  frames  of  the  horizontally  and  randomly  polarized 
images  in  a  single  data  file  is  roughly  8  seconds.  Within  this  time  scale,  the  temperature 
distribution  of  the  ship  surfaces  has  not  changed  significantly.  If  we  translate  the  first 
image  frame  of  each  polarization  in  a  data  file  into  the  equivalent  blackbody  temperature 
using  the  calibration  curves,  we  should  expect  the  blackbody  temperature  arrays  for  the 
three  different  polarizations  to  be  almost  similar,  that  is,  a  pixel  in  one  polarization  frame 
should  have  the  same  equivalent  blackbody  temperature  as  the  corresponding  pixel  in 
another  polarization  frame.  Hence,  if  we  take  the  difference  between  the  blackbody 
temperature  arrays  of  any  two  polarizations,  the  elements  of  the  resultant  temperature 
difference  matrix  would  be  mostly  zeros.  In  another  way,  the  average  value  of  all  the 
elements  in  the  resultant  temperature  difference  matrix  should  be  nearly  zero.  Table  4.1 
below  shows  the  average  values  of  the  elements  of  the  resultant  temperature  difference 
matrices  when  the  blackbody  temperature  arrays  of  the  three  different  polarizations  of  a 
single  data  file  are  compared  to  each  other.  Column  one  gives  the  average  temperature 
difference  between  horizontal  and  vertical  polarizations;  column  two  ,  horizontal  and 
random  polarizations;  and  column  three  vertical  and  random  polarizations. 
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Table  4.1  Mean  pixel  temperature  difference  between  polarization  cases  of  an  image 

file 


Filename 

Average  temperature  difference  (K) 

H-V 

H-U 

V-U 

Base09.pte 

-0.053289222 

-6.0796106 

-6.0263214 

BaselO.pte 

-0.063456130 

-4.8232400 

-4.7597838 

Basel  l.pte 

0.094446067 

-5.2094861 

-5.3039322 

BaselS.pte 

0.029233603 

-4.7203123 

-4.7495459 

Basel6.pte 

0.23777563 

-4.1678803 

-4.4056560 

Base28.pte 

0.55287224 

-1.4238610 

^  -1.9767332 

Bases  O.pte 

0.21753705 

-1.3070561 

-1.5245932 

Base32.pte 

-0.45393654 

-3.3308119 

-2.8768753 

We  can  see  that  the  temperature  differences  between  corresponding  pixels  of  the 
polarized  image  and  unpolarized  image  are  significant.  If  this  difference  were  solely  due 
to  the  contrast  improvement  of  the  polarizing  filters,  we  should  expect  opposite  signs  for 
the  second  and  third  columns.  Clearly,  this  is  not  the  case  and  it  is  likely  that  calibration 
drifts  are  present.  Since  we  do  not  know  exactly  what  the  calibration  drifts  are,  a  relative 
calibration  is  required  in  order  for  us  to  relate  the  measured  ship  skin  temperatures  to  the 
recorded  thermal  levels.  The  procedure  for  a  relative  measurement  is  described  in  section 
10  of  the  AGA  manual  [Ref.  10  ].  The  relevant  extract  is  reproduced  in  Appendix  O. 
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To  perform  the  relative  calibration,  we  use  the  following  formula  from  the  AGA 
manual : 


where  Iq  is  the  relative  calibrated  thermal  level, 

If  is  the  calibrated  thermal  level  of  a  reference  point  on  the  object  whose 
temperature  is  known, 

Alor  is  the  thermal  level  difference  from  the  reference  level  (T), 

To  is  the  atmospheric  transmittance  between  object  and  scanner  and 
So  is  the  object  emissivity. 

Since  the  values  of  the  elements  of  the  matrix  that  we  derived  are  radiance  values 
received  at  the  scanner  aperture,  these  would  have  to  be  converted  into  equivalent 
blackbody  temperatures.  The  routine  RADTOT.PRO  converts  radiance  to  equivalent 
blackbody  temperature  by  linear  interpolation  in  three  iterations.  Since  the  radiation- 
temperature  relation  is  actually  non-linear,  a  good  initial  guess  of  the  temperature  range  is 
required  to  achieve  good  conversion  accuracy.  For  the  temperature  range  of  the  ship,  the 
conversion  accuracy  is  0.001  W/m^. 

If  we  now  use  the  converted  blackbody  temperature  of  the  pixels  for  the 
recalibration,  then  to  and  So  are  simply  equal  to  one  in  equation  4.3  since  these  are  the 
equivalent  blackbody  temperatures  at  the  aperture  of  the  scanner.  In  addition,  any  pixel 
can  be  used  as  the  reference  point  since  now  they  each  have  a  known  equivalent 
blackbody  temperature. 

4.  Unit  of  Comparison 

The  comparison  can  be  done  using  radiance,  temperature  or  thermal  level 
differences.  However,  thermal  level  differences  do  not  give  a  good  indication  when  we 
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compare  across  polarization  cases.  For  example,  a  thermal  level  difference  of  1  in  the 
horizontal  polarization  case  may  be  due  to  0.3  W/m^  difference  in  radiance  while  a 
difference  of  3  isotherm  units  in  the  unpolarized  case  is  just  0. 1 5  W/ml  Hence,  doing  the 
comparison  in  temperature  units  is  the  most  convenient  since  the  recalibrated  thermal 
levels  of  the  image  pixels  are  easily  translated  to  temperature. 

B.  DISCUSSION  ON  RESULTS 

1.  Statistical  Measures 

In  an  ideal  case,  the  temperature  difference  between  each  pixel  of  the  model  and 
the  recorded  image  would  be  zero.  More  realistically  speaking,  we  would  expect  a 
distribution  of  pixels  with  a  mean  of  close  to  zero  and  a  small  standard  deviation.  Figure 
4.1  below  shows  the  pixel  distributions  for  a  sample  file.  From  the  figure,  a  large  number 
of  pixels  are  clustered  close  to  the  zero  temperature  difference  mark.  However,  there  are 
a  small  percentage  of  pixels  that  register  a  big  temperature  difference.  The  origins  of 
these  outlying  pixels  are  likely  due  to  two  limitations. 

The  first  limitation  is  inherent  in  the  geometric  representation  of  the  ship.  As  our 
geometric  model  is  a  simplified  representation  of  the  ship,  not  all  calculated  pixels  will 
have  a  corresponding  pixel  in  the  image  or  vice  versa.  Secondly,  the  alignment  of  the 
mask  is  highly  sensitive,  that  is,  down  to  pixel  resolution.  At  this  level  of  resolution,  it  is 
difficult  to  achieve  a  'correct'  one  to  one  correspondence  in  the  ship  pixels.  The  impact  of 
all  this  is  that  some  pixels  will  be  compared  with  the  backgroimd  pixels  resulting  in  a  big 
temperature  difference.  It  has  the  effect  of  driving  up  the  mean  and  standard  deviation 
especially  when  the  number  of  ship  pixels  is  small. 

In  order  to  obtain  a  better  characterization  of  the  pixel  distribution,  two  additional 
statistical  measures  were  calculated.  These  are  the  modal  temperature  difference  and  the 
fraction  of  pixels  falling  within  a  specified  range  firom  the  modal  temperature  difference. 
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The  modal  temperature  difference  has  the  advantage  that  it  is  not  affected  by  the  outlying 
pixels.  However,  it  is  more  susceptible  to  freak  occurrence  of  peaks  in  a  distribution. 
Therefore,  both  sets  of  statistical  measures  are  needed  to  gain  a  fair  interpretation  of  the 
results. 


Figure  4.1  Pixel  distribution  of  difference  in  equivalent  blackbody  temperature 
between  model  and  recorded  image  for  file  Basel O.pte.  The  small  percentage  of 
pixels  at  the  outlying  areas  tend  to  drive  up  the  mean  and  standard  deviation 
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2.  Results  of  Comparison 


A  total  of  37  data  files  were  analyzed  and  compared.  The  results  of  each 
comparison  are  tabulated  in  Table  4.2.  Most  of  the  cases  analyzed  showed  a  modal 
temperature  difference  of  less  than  IK  between  model  and  image.  The  average 
temperature  differences  for  these  cases  were  also  less  than  IK.  The  modal  temperature 
difference,  in  general,  is  lower  than  the  average  temperature  difference. 

There  are  5  files  where  the  modal  temperature  difference  is  much  higher  than  the 
average  temperature  difference.  The  modal  temperature  differences  for  the  3  polarization 
cases  in  these  files  were  very  different  as  well.  In  contrast,  the  average  temperature 
differences  for  the  3  polarization  cases  were  quite  consistent.  It  is  likely  that  there  are 
occurrences  of  freak  peaks  in  the  distribution  making  the  modal  temperature  difference 
unreliable. 

There  is  one  file  that  has  average  and  modal  temperature  differences  of  more  than 
2K  for  all  3  polarization  cases.  This  is  file  BaseOT.pte.  The  total  number  of  pixels  for 
each  frame  of  this  file  is  32.  It  is  likely  that  there  are  too  few  good  pixels  to  give  a 
reasonable  comparison.  Withthe  exception  of  file  BaseOT.pte,  there  is  moderate 
agreement  between  the  model  and  recorded  images  with  45-50%  of  the  pixels  falling 
within  1 .75K  of  the  measured  apparent  temperature. 
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Table  4.2  Comparison  of  temperature  between  model  and  recorded  images 


Measure  of  spread 

Standard 

deviation 

00 

00 

00 

3.033 

o 

00 

3.007 

9Z0‘Z 

3.617 

3.578 

3.730 

4.894 

3.552 

2.906 

Fraction  of  pixels 
within  0.75  K  of 
modal  temp, 
difference 

0.437 

0.459 

0.602 

0.506 

0.526 

0.553 

0.476 

OO 

OO 

d 

0.349 

0.213 

0.434 

0.533 

Fraction  of  pixels 
within  0.25  K  of 
modal  temp, 
difference 

0.266 

ZZZ'O 

L9Z0 

0.219 

OO 

OO 

d 

0.288 

0.280 

0.293 

0.171 

0.128 

0.213 

0.277 

Measure  of  centre  of 
distribution 

Mean 

temperature 

difference 

2.391 

1.223 

_ 

-0.282 

1.332 

0.218 

-0.569 

1.136 

OO 

1.388 

0.212 

0.136 

-0.511 

Modal 

temperature 

difference 

oo 

1.182 

0000 

1.566 

6600 

-0.476 

-0.389 

69  ro 

0.016 

-2.726 

-0.274 

-0.787 

Polarisation 

case 

> 

X 

> 

D 

X 

> 

D 

X 

> 

D 

Number 
of  pixels 

00 

yn 

O 

CN 

OO 

File  no. 

(base*.pte) 

o 

(N 

O 

o 

o 
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Modal  Mean  Fraction  of  pixels  Fraction  of  pixels  Standard 

temperature  temperature  ivithin  0.25  K  of  within  0.75  K  of  deviation 
difference  difference  modal  temp.  modal  temp. 

difference  difference 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


The  simplified  surface  radiosity  equation  is  the  basis  for  the  modeling.  It  states 
that  for  each  pixel  the  radiance  from  the  surface  is  the  sum  of  the  emitted  radiance  and 
reflected  radiance.  Since  we  have  made  the  Lambertian  and  graybody  approximations, 
the  emitted  radiance  accounts  for  almost  95%  of  the  radiance  before  atmospheric 
attenuation.  From  the  comparison  of  apparent  temperatures  between  the  model  and  the 
images,  55%  of  the  pixels  have  temperature  deviations  of  more  than  1 .75K  from  the 
images.  For  the  range  of  temperatures  read  by  the  thermistors,  this  deviation  in  radiance 
is  of  the  order  of  about  3  W/m^-sr.  We  can  see  from  the  SeaRad  outputs  (Appendix  K) 
that  the  average  sea  and  sky  radiance  incident  on  the  pixel  surface  is  of  the  order  of  30 
wW-sr  and  10  wW-sr  respectively.  This  gives  the  reflected  radiance  to  be  in  the  order 
of  1  W/m^-sr.  Therefore,  the  origin  of  the  1 .75K  temperature  deviations  for  55%  of  the 
pixels  cannot  to  be  due  to  the  reflected  radiance.  A  second  possibility  for  the  deviation  is 
the  path  radiance  due  to  path  length  error.  The  path  radiance  is  of  the  order  of  0.01  W/m  - 
sr  per  meter  increase  in  range.  It  would  require  a  range  error  of  more  than  300m  to 
accoxmt  for  the  3  W/m^-sr  deviation.  Again,  this  is  unlikely  to  be  the  source  of  deviations 
since  the  range  or  elevation  error  is  not  expected  to  be  more  than  10m.  Therefore,  we  can 
conclude  that  the  mostly  likely  source  of  the  error  is  the  surface  temperature  distribution 
used  in  the  model.  It  does  not  sufficiently  represent  the  actual  ship  skin  temperature 
distribution. 

One  of  the  problems  with  modeling  the  temperature  distribution  of  the  ship  is  that 
there  are  too  few  thermistor  readings.  As  a  result,  each  plane  is  assumed  to  have  a 
constant  temperature  associated  with  the  nearest  thermistor.  Given  the  limitation  on  the 
number  of  thermistors,  one  option  is  to  perform  sensitivity  analysis  with  the  model.  This 
will  help  in  a  better  imderstanding  of  the  actual  temperature  distribution  of  the  ship.  One 
approach  is  to  perform  a  plane  by  plane  analysis.  Additional  planes  can  be  added  to  give 
a  better  distribution  where  necessary.  Another  approach  may  be  to  consider  one  or  two 
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dimensional  heat  flow  for  each  plane  so  as  to  arrive  at  a  temperature  distribution  for  the 
pixel  surface  elements. 

As  discussed  earlier,  the  main  contributions  to  the  observed  apparent  radiance  are 
the  emitted  radiance  from  the  ship  surfaces  and  the  path  radiance.  Since  we  have  assumed 
these  two  contributions  to  be  randomly  polarized,  the  apparent  radiance  that  we  derived 
from  our  ship  model  would  show  no  significant  polarization  features.  When  our  model  is 
compared  to  the  measured  images  for  the  three  polarization  cases,  the  match  is 
comparable  in  all  three  polarizations.  The  average  overall  modal  temperature  differences 
are  0.078K,  -0.093K  and  -0.093K  for  horizontal,  vertical  and  random  polarization 
respectively.  The  spreads  of  the  distributions  for  all  three  polarizations  are  also 
comparable  at  40-45%  of  the  pixel  within  0.75K  of  the  mode.  This  shows  that  the 
assumptions  used  in  our  model  are  justified. 

For  the  whole  modeling  process,  the  running  of  SeaRad  code  is  the  most  time 
consuming.  This  is  because  each  input  file  is  copied  to  the  standard  input  file  to  run  the 
SeaRad  code  and  the  standard  output  file  is  then  copied  to  an  output  file.  As  each  pixel 
requires  one  input  file,  it  takes  approximately  one  hour  to  run  a  frame  consisting  of  360 
pixels.  It  also  takes  up  a  lot  of  memory  space  to  store  the  input  and  output  files.  The 
SeaRad  input  data  requirement  can  be  broken  down  into  two  types  of  data.  The  first  type 
is  the  meteorological  data.  This  data  does  not  change  for  one  frame  and  is  the  same  for  all 
input  files  of  that  frame.  The  second  type  is  the  geometry  input.  This  is  the  part  that 
changes  for  each  pixel.  If  it  were  possible  to  modify  the  SeaRad  code  to  accept  multiple 
geometry  input,  then  this  would  significantly  reduce  the  processing  time  as  well  as  the 
time  spent  copying  files.  It  will  make  the  SeaRad  code  a  much  more  practical  tool. 


42 


APPENDIX  A.  DIMENSIONS  OF  PLANES  IN  GEOMETRIC 
MODEL  AND  GENERAL  DATA  ON  RA^  POINT  SUR 


A.  DIMENSIONS  OF  PLANES  IN  GEOMETRIC  MODEL 


The  geometric  representation  of  the  ship  is  made  up  of  23  vertical  planes.  An 
illustration  of  the  ship  model  is  at  Figure  2.2  in  Chapter  2.  The  table  below  summarises 
the  dimensions  of  each  plane  in  the  model. 


Table  Al.  Summary  Of  Planes  In  Ship  Model 


Plane 

Length 

Height 

Matched  to 

Description 

No. 

(m) 

(m) 

thermistor 

1 

12.6 

4.5 

10 

Port  hull  mid  section 

2 

6.9 

4.5 

4 

Port  hull  below  stack 

3 

7.5 

2.7 

2 

Port  hull  aft  section 

4 

9.6 

3.3 

8 

01  deck  port  side 

5 

3.6 

3.6 

14 

Stack  port  side 

6 

14.4 

2.7 

7 

Starboard  hull  aft  section 

7 

12.6 

4.5 

16 

Starboard  hull  mid  section 

8 

9.6 

3.3 

15 

01  deck  starboard  side 

9 

3.3 

8 

02  deck  portside 

00 

3.3 

15 

02  deck  starboard  side 

11 

6.9 

1.8 

4 

Electrical  switchboard  room 

43 


Plane 

No. 


12 

13 

14 


15 

16 
T? 
Ts 

19 

iT 


Length 

(m) 


Height 

(m) 


3.6 


3.6 


9.6 


2.7 


4.2 


1.8 


5.4 


1.8 


9.6 


3.3 


7.2 


3.3 


1.8 


9.6 


7.2 


1.8 


12.9 


3.6 


3.3 


3.3 


3.6 


5.4 


Matched  to 
thermistor 


Description 


12 


2 


4 


7 


6 


6 


5 


15 


15 


14 


11 


starboard  side 
Stack  starboard  side 
Stem 

Electrical  switchboard  room 
aft 

00  deck  aft 
01  deck  aft 
02  deck  aft 
Stack  aft 
01  deck  fore 
02  deck  fore 
Stack  fore 
Port  bow 


23 


12.9 


5.4 


16 


Starboard  bow 


APPENDIX  B.  THERMISTOR  LOCATIONS  ON  RA^  POINT  SUR 


The  thermistor  location  is  described  at  Table  Bl.  Figure  B1  shows  the  locations. 


Table  Bl.  Location  of  thermistors 


Channel  No. 

Location 

1 

Reference  25 

2  • 

Aft  port 

3 

Aft  water  probe 

4 

Aft  of  stack 

.  5 

Aft  on  stack 

6 

Aft  boathouse 

7 

Aft  starboard 

8 

Port  boathouse 

9 

Reference  25 

10 

Port  below  stack 

11 

Port  bow 

12 

Starboard  stack 

13 

Air 

14 

Port  stack 

15 

Bow  boathouse 

16 

Starboard  bow 
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®  d)  CD 


Figure  Bl.  Location  of  thermistors 
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APPENDIX  C.  AVERAGE  SEA  AND  SKY  RADIANCE 
APPROXIMATIONS 


A.  AVERAGE  SKY  RADIANCE  PATH 

The  sky  radiance  may  be  considered  to  be  the  sum  of  the  sky  background 
emission  and  the  path  radiance  for  a  cloudless  sky.  It  is  reasonable  to  assume  that  the  sky 
background  is  isotropic  and  therefore,  sky  background  emission  is  independent  of  look 
angle.  However,  path  radiance  would  increase  with  increasing  zenith  angle  towards  the 
horizon  since  the  path  length  is  longest  at  the  horizon.  Hence,  the  path  representing  the 
average  sky  radiance  is  the  same  path  as  that  for  the  average  path  radiance. 

To  obtain  the  average  path  radiance,  we  consider  the  geometry  of  the  atmospheric 
layer  and  the  earth  surface  in  Figure  Cl  below.  SeaRad  assumes  an  atmospheric  layer  of 
100  km  from  the  earth  surface. 


Figure  Cl.  Relationship  between  zenith  angle,  6  and  angles  a  and  P 
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The  angles  0,  a  and  p  are  simply  related  by  the  sine  rule  : 


sin(180-^)  _  siny^  _  sina 
i?  +  100  ~  R  P 


(C.1) 


Therefore,  we  find  that  the  path  length,  P  is  a  function  of  0.  The  average  path 
length  is  then  given  by  : 


y  nil 

{P)  =  -  jP(0)dd 

^  n 


(C2) 


Performing  the  integral  numerically,  we  obtain  an  average  path  length  of  208  km, 
which  corresponds  to  a  slant  path  to  space  at  a  zenith  angle  of  65.3  degrees. 


B.  AVERAGE  SEA  RADIANCE  PATH 

Assuming  the  sea  to  be  a  flat  calm  surface,  the  sea  radiance  may  be  simplified  as  : 

Lsea  “  Bsea  (T)-p(0i)[Bsea(T)-Lsky]  (C.3) 


where  p  is  the  reflectance  given  by  the  Fresnel  equations , 

0i  is  the  angle  of  incidence  or  elevation, 

Bsea  is  the  blackbody  emission  of  the  sea  surface  at  temperature  T  and 
we  have  assumed  here  that  Lsky  is  the  average  sky  radiance. 

We  can  see  that  Lsea  is  also  a  function  of  the  zenith  angle.  Since  Bsea  and  Lsky  are  both 
constant,  the  path  representing  the  average  sea  radiance  is  the  same  path  as  that  for  the 
average  reflectance. 
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From  the  Fresnel  equations,  the  E-vector  amplitude  reflection  ratios  are  : 


cosdj  -  n,  cosO, 
n,cos^,  +n,cos^, 


(C.4a) 


n,  cosOj  -rij  cosO,  (C.4b) 

HjCosOi +n,cos6j 


Taking  nj  =1  and  nt  =1 .33  and  making  use  of  the  Snell's  law,  we  can  find  the  reflectance, 
p  as  a  function  of  6i ,  where  : 


p(0i )  =  ri^  +  r/ 


(C.5) 


Hence,  the  average  reflectance  is  : 

^  0 

Again,  performing  the  integral  numerically,  we  found  the  average  reflectance  of  the  water 
surface  to  be  0.208.  This  value  corresponds  to  the  slant  path  to  sea  with  a  zenith  angle  of 
118.1  degree. 
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APPENDIX  D.  LIST  OF  PROGRAMS 


Part  1 


Start.pro 

procedure 

Main  procedure  to  extract  meteorological  data  and 
determine  aspect  angle  and  pixel  size 

Readtime 

function 

Reads  the  local  time  (PST)  from  *.pte  files 

Readmet.pro 

function 

Reads  meteorological  data  and  returns  a  met 
vector  containing  all  the  required  information 

Leiajpte.pro 

function 

Extracts  the  images,  information  or  time  vector 

Planes.pro 

ftmction 

Selects  the  planes  that  are  visible  for  a  given 
aspect  angle 

Planein.pro 

procedure 

Calculates  the  geometry  inputs  for  SeaRad  files  for 
each  pixels 

Seain.pro 

procedure 

Writes  SeaRad  input  files  between  two  altitude 

Avein.pro 

procedure 

Writes  SeaRad  input  files  for  average  sea  and  sky 
paths 

Posn.pro 

function 

Calculates  the  position  of  each  node 

Rposn.pro 

function 

Calculate  the  position  of  each  node  after  rotation 
to  a  given  aspect  angle 

Vector  .pro 

function 

Calculates  the  vector,  V,  from  the  sensor  to  the 
node  position 

Vsea.pro 

function 

Calculates  the  vector  from  the  nodes  to  tiie  sea 

Vmag.pro 

function 

Calculates  the  magnitude  of  the  V  vector 

Zenith.pro 

function 

Calculates  the  zenith  angle  at  the  sensor  for  each 
node 

Part2 

Seal 

DOS 
batch  file 

Copy  all  input  files  in  Case  1 /Inputs  directory  to 
tape  5  and  call  XI 

XI 

DOS 
batch  file 

Execute  SeaRad  and  copy  output  file  to 

Case  1 /Outputs  directory 
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File  Name 

Type 

Function 

Parts 

Continue.pro 

procedure 

Main  procedure  to  continue  the  analysis  after 
running  SeaRad 

Readtemp.pro 

function 

Reads  the  temperature  matrix 

Seaout.pro 

function 

Reads  the  output  files  from  SeaRad  and  returns  a 
matrix  with  the  transmittance  and  path  radiance 

Aveout.pro 

function 

Reads  the  output  files  from  SeaRad  and  returns  a 
vector  of  average  sky  and  sea  radiance 

Emission.pro 

function 

Calculates  the  in-band  radiance 

Place  .pro 

function 

Aligns  all  the  radiance  values  of  planes  into  one 
matrix 

Block.pro 

procedure 

Determines  areas  of  overlap  between  planes 

Overlap.pro 

function 

Corrects  the  radiance  values  in  the  overlapped  area 

Radtot.pro 

function 

Calculates  the  equivalent  blackbody  temperature 
from  radiance 

Convtemp.pro 

function 

Converts  radiance  matrix  to  temperature 

Ttotl.pro 

function 

Calculates  thermal  level  from  temperature  based 
on  calibration  constants 

Dltotl.pro 

function 

Converts  digital  level  to  thermal  level 
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APPENDIX  E.  STRUCTURE  OF  CONVERTED  DATA  FILE 


The  converted  files  have  the  following  properties  : 

-  all  polarisation  cases  for  one  test  case  are  combined  into  one  file  rather  than 
having  two  separate  files 

-each  image  is  a  125  x  256  integer  array 

-  all  essential  information  is  put  into  a  single  integer  matrix 

The  structure  of  the  file  and  the  format  for  header  information  is  tabulated  below. 


Table  El.  Structure  Of  A  Basic  (.Pte)  File 


Variable 

Type 

Name  of  the  original  (.PTW)  file  containing 
horizontal  and  vertical  polarisation  cases 

string(20) 

Information  vector 

integer  (1,13) 

Frame  time  :  horizontal  polarisation 

integer  (4,  #  frames  hor.case) 

Frame  time  :  vertical  polarisation 

integer  (4,#  frames  ver.case) 

Frame  time  :  unpolarised  polarisation 

integer  (4,#  frames  impol.case) 

Image  set  (frames)  :  horizontal  polarisation 

integer  (125,256,#  frames  hor.case) 

Image  set  (frames)  :  vertical  polarisation 

integer  (125,256,#  frames  ver.case) 

Image  set  (frames) :  unpolarised  polarisation 

integer  (125,256,#  frames  impol.case) 

Table  E2.  Information  Vector 


Line 

Data 

Example 

0 

Day 

9 

1 

Month 

4 

2 

Year 

96 

3 

Hour  (first  frame) 

12 

55 


Line 

Data 

Example 

4 

Minute  (first  frame) 

45 

5 

Seconds  (first  frame) 

27 

6 

Himdreds  (first  frame) 

77 

7 

Number  of  polarisation  cases 

o 

J 

8 

Number  of  pixels  per  line 

125 

9 

Number  of  lines 

256 

10 

Number  of  frames  :  horizontal  polarisation 

27 

11 

Number  of  frames  :  horizontal  polarisation 

25 

12 

Number  of  frames  :  horizontal  polarisation 

29 

56 


36tn  £c..Mav2  98  ■ 

-I 
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APPENDIX  G.  STORAGE  FORMAT  FOR  METOC  DATA 


There  were  two  meteorological  towers  installed  onboard  RA^  Point  Sur  for 
recording  the  meteorological  data.  The  storage  format  for  the  two  towers  are  slightly 
different  as  the  METOC  2  system  has  two  additional  temperature  and  relative  humidity 
sensors.  The  two  formats  shall  be  described  below. 

A.  METOC  1  DATA 

These  data  files  are  identified  by  the  filename  "eop*.txt"  where  *  is  the  GMT  day 
of  the  year.  Each  file  is  started  at  00:00  GMT  and  the  data  are  averaged  every  30  seconds 
and  recorded.  For  each  line  of  a  record,  the  various  readings  from  the  sensors  are 
separated  by  spaces  across  a  line  of  text.  At  the  end  of  the  record,  a  carriage  return  and 
line  feed  character  starts  the  next  record  on  a  new  line.  The  order  of  the  sensor  readings  is 
as  follows : 

0 1  number  of  readings  in  this  average,  nn 
02  year,  yy 
03  Julian  date,  ddd 
04  hours  &minutes,  hhmm 
05  seconds,  ss 

06  relative  wind  speed  (m/s),  mmm 

07  relative  wind  direction  (degrees),  ddd 

08  ship  speed  (knots),  kk.k  -  information  not  available 

09  ship  direction  (degrees),  ddd.d  -  information  not  available 

10  true  wind  speed  (m/s),  mmm 

1 1  true  wind  direction  (degrees),  ddd 

12  Air  temperature  (degrees  C),  tt.t 

13  relative  humidity  (%),  pp 

14  pressure  (millibars),  pppp.p 
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15  sea  surface  temperature  (degrees  C),  tt.t 

16  GPS  time  (GMT),  hhmmss 

17  GPS  latitude  (degrees  &  minutes),  ddmm.m 

18  GPS  north/south  indicator,  either  "N'  or  'S' 

19  GPS  longitude  (degrees  &  minutes),  ,  ddmm.m 

20  GPS  east/west  indicator,  either  'E'  or  'W 

21  GPS  speed  over  ground  (knots),  kk.k  -  information  not  available 

22  GPS  course  over  ground  (degrees),  ddd.d  -  information  not  available 

23  GPS  antenna  height  (m),  mm.m 

A  sample  of  the  data  is  below : 

30  1996  99  2352  32  2  138  0  090  5  288  16.5  75  1016.9 

16.8  000035  3311.2  N  11728.2  W  6.3  124.3  46.3 

B.  METOC  2  DATA 

METOC  2  data  files  are  identified  by  the  filename  "eop30sec.*  "  where  *  is  the 
GMT  day  of  the  year.  The  recording  format  is  the  same  as  METOC  1  files  except  that 
there  are  4  additional  sensor  readings  : 

01  number  of  readings  in  this  average,  nn 
02  year,  yy 
03  Julian  date,  ddd 
04  hours  &minutes,  hhmm 
05  seconds,  ss 

06  relative  wind  speed  (m/s),  mmm 

07  relative  wind  direction  (degrees),  ddd 

08  ship  speed  (knots),  kk.k  -  information  not  available 

09  ship  direction  (degrees),  ddd.d  -  information  not  available 

10  true  wind  speed  (m/s),  mmm 

1 1  true  wind  direction  (degrees),  ddd 


60 


12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 

27 

28 


Air  temperature  (degrees  C),  tt.t 

relative  humidity  (%),  pp 

pressure  (millibars),  pppp.p 

sea  surface  temperature  (degrees  C),  tt.t 

Air  temperature  A  (degrees  C),  tt.t 

relative  humidity  A  (%),  pp 

Air  temperature  B  (degrees  C),  tt.t 

relative  humidity  B  (%),  pp 

"***",  Campbell-GPS  data  delimiter 

GPS  time  (GMT),  hhmmss 

GPS  latitude  (degrees  &  minutes),  ddmm.m 

GPS  north/south  indicator,  either  'N'  or  'S' 

GPS  longitude  (degrees  &  minutes),  ,  ddmm.m 
GPS  east/west  indicator,  either  'E'  or  'W 

GPS  speed  over  groxind  (knots),  kk.k  -  information  not  available 
GPS  course  over  ground  (degrees),  ddd.d  -  information  not  available 
GPS  antenna  height  (m),  mm.m 


A  sample  is  below  : 

30  1996  99  2352  32  2  138  0  090  5  288  16.5  75  1016.9 

16.8  0.0  0.0  0.0  0.0***  000035  331 1.2  N  11728.2  W  6.3  124.3 

46.3 
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APPENDIX  H.  EXTRACTED  INFORMATION  OF  ANALYSED  FILES 


Line  of 
Sight 
(LOS) 
(deg) 

246.2 

246.2 

246.2 

241.9 

241.9 

241.9 

243.2 

245.0 

245.0 

245.0 

247.9 

256.4 

247.9 

Wind 
azimuth 
relative  to 
LOS  (deg) 

37.8 

37.8 

32.8 

49.8 

63.1 

49.2 

57.0 

57.0 

51.0 

36.1 

27.6 

26.1 

Average 

wind 

speed 

(m/s) 

5.2 

5.20 

5.20 

5.20 

5.20 

5.20 

5.20 

5.20 

1 

5.20 

Current 

wind 

speed 

(m/s) 

o 

o 

O 

O 

O 

O 

to 

009 

009 

5.67 

O 

O 

o 

o 

O 

O 

O 

o 

1 

O 

O 

to 

AMP 

cn 

cn 

m 

VO 

Sea  temp. 

(“o 

0891 

16.80 

17.20 

15.87 

o 

oo 

14.80 

15.07 

15.10 

15.10 

o 

00 

14.20 

14.20 

15.30 

Range 

(km) 

■ 

1.31 

1.31 

1.31 

2.07 

2.07 

2.07 

2.07 

2.63 

2.63 

2.63 

3.56 

«  WD 

i»  s  ^ 

«  w 

57.2 

80.2 

163.2 

29.7 

71.5 

164.8 

9.9 

32.9 

54.9 

00 

O 

o 

10.9 

GMT 

(hrs) 

m 

CN 

o 

CN 

o 

CN 

O 

o 

CN 

o^ 

o 

CN 

r- 

o 

CN 

20.79 

20.86 

21.07 

21.13 

21.18 

21.25 

21.33 

21.47 

Day  of 
year 

o 

o 

O 

o 

o 

o 

o 

o 

t-H 

100 

o 

o 

100 

O 

o 

O 

O 

O 

O 

100 

File  no. 
(base*.pte) 

o 

CN 

o 

cn 

o 

o 

05 

90 

07 

08 

60 

O 

- 

CN 

cn 
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21.52  31.9  3.56  15.30  6  5.00  5.20  26.1  247.9 
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File  no.  Day  of  GMT  Aspect  Range  Sea  temp.  AMP  Current  Average  Wind  Line  of 

(base*.pte)  year  (hrs)  angle  (km)  ("C)  wind  wind  azimuth  Sight 

(deg)  speed  speed  relative  to  (LOS) 

(m/s)  (m/s)  LOS  (deg)  (deg) 
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APPENDIX  I.  NAVY  AEROSOL  MODEL  AND  AIR  MASS 

PARAMETER 


A.  NAVY  AEROSOL  MODEL  [Ref.  5] 


The  US  Navy  Aerosol  Model  (NAM)  was  developed  to  compute  the  aerosol 
extinction  coefficients  for  visual  and  infrared  wavelengths  over  a  marine  boundary 
environment.  It  is  built  from  meteorological  and  aerosol  size  spectra  databases  obtained 
near  the  surface  layer  of  several  oceanic  environments. 


The  aerosol  size  distribution  of  the  NAM  is  attributed  to  3  distinct  aerosol 
production  mechanisms.  The  first  mechanism  is  air  mass  driven  while  the  second  and 
third  mechanisms  are  wind  driven.  Correspondingly,  the  NAM  aerosol  size  distribution  is 
a  sum  of  3  independent  log-normal  distribution  functions.  This  is  given  as  follows  ; 


,  .  dN  ^  .  I 
n{r)  =  —  =  ^-A,QXid 

dr  M  /  ' 


V  1-^/. 


(1.1) 


where  n(r)  is  the  number  of  particles  per  cm^  per  radius  r,  in  microns  (pm), 
k  is  the  visibility  scaling  factor, 

ri  are  the  modal  radii  (  0.03pm,  0.24  pm  and  2.0pm)  measured  at  80%  RH,  and 
A 1  =  2000  ,  where  P  is  the  air  mass  parameter, 

A2  =  5.866(Vave-2.2)  ,  where  Vave  is  the  24-hour  average  wind  velocity. 

As  =  10  ,  where  V  is  the  instantaneous  wind  velocity, 

/  =  [(2-Rh/100)/6(l-Rh/100)f^ 


B.  AIR  MASS  PARAMETER  [Ref.  5] 


Of  the  3  amplitudes  Ai  to  A3  given  above,  the  largest  and  the  most  difficult  to 
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measure  is  the  first  one,  that  is,  the  air  mass  driven  component.  The  air  mass  parameter,  P 
is  an  integer  ranging  from  1  for  clean  oceanic  air  to  10  for  polluted  continental  or  coastal 
air. 


There  are  several  methods  to  obtain  the  air  mass  parameter.  These  methods  range 
from  measuring  the  atmospheric  Radon-222  gas  concentration  to  computing  air  mass 
trajectories  to  remote  sensing  using  satellite.  Each  method  has  its  advantages  and 
disadvantages. 

Traditionally,  the  Radon-222  gas  concentration  method  is  the  most  reliable. 
Radon-222  is  an  inert  radioactive  gas  with  a  half-life  of  3.82  days  produced  by  the  decay 
of  Radium-236.  Radium-236  is  found  in  various  quantities  in  the  soils  and  rocks  of  the 
earth's  crust.  After  formation  in  the  groimd,  the  radon  diffuses  into  the  atmosphere  and 
becomes  attached  to  aerosol  particles  that  are  transported  by  the  prevailing  wind  over  the 
ocean.  Since  radon  is  mainly  produced  on  land,  it  is  used  as  an  air  mass  tracer  to 
characterise  the  air  mass  history.  The  air  mass  parameter,  P  is  then  simply  given  by  the 
equation : 


(1.2) 

where  INT  truncates  to  the  integer  value  and  Rn  is  the  atmospheric  Radon-222  gas 
concentration  in  picocurries  per  cm^.  Strictly  speaking,  the  concentration  of  the  radon  in 
different  areas  will  not  be  the  same.  This  is  because  it  depends  on  the  soil  conditions, 
atmospheric  conditions,  radium  deposit  in  the  crust  and  strength  of  the  radium  source. 
Therefore,  it  is  the  relative  change  in  radon  gas  concentration  that  is  meaningful. 
Equation  1.2,  however,  uses  the  absolute  radon  concentration.  This  is  because  the  factor 
of  4  was  established  for  the  region  at  San  Nicolas  Island.  This  factor  will  have  to  be 
modified  for  other  geographical  regions. 


P  =  INT\ 


Rn 

~4~ 
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APPENDIX  J.  SEARAD  INPUT 


SeaRad  input  requires  the  submission  of  a  series  of  a  minimum  of  5  cards.  The 
number  of  cards  required  depends  on  the  mode  of  execution.  There  are  a  total  of  19  cards 
of  which  5  are  the  main  cards  and  the  rest  are  optional  cards.  A  complete  description  of 
each  card  can  be  found  in  the  file  "ModtranlO"  that  comes  with  program.  However,  this 
does  not  include  the  modifications  done  in  SeaRad.  The  file  "notes.txt"  logs  all  the 
modifications  and  should  be  read  in  conjunction  with  the  original  card  description.  This 
appendix  will  only  explain  the  parameters  that  are  required  in  the  cards  needed  in  this 
thesis. 


CARD  1  (line  one)  -  Format  (1315,  F8.3,  F7.2) 

MODEL,  ITYPE,  lEMSCT,  IMULT,  Ml,  M2,  M3,  M4,  M5,  M6,  MDEF,  IM,  NOPRT, 
TBOUND,  SALE 

F  -  use  LOWTRAN  7 

MODEL  -  selects  one  of  the  6  geographical  model  atmospheres  or  specifies  that 

user-defined  meteorological  data  are  to  be  used 
:  mid-latitude  summer  (2)  was  used 
ITYPE  -  indicates  the  type  of  atmospheric  path 

.•  slant  path  to  space  (3)  was  used  for  average  sea  and  sky  radiance 
calculations  and  slant  path  between  two  altitudes  (2)  was  used  for  pixel¬ 
sensor  atmospheric  compensation 
lEMSCT  -  selects  the  execution  mode  of  the  program 

.•  radiance  mode  with  solar  scattering  (2)  used 
IMULT  -  determines  execution  with  multiple  scattering 
.•  multiple  scattering  (1)  was  used 
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Ml  ,M2,M3,M4,M5,M6  -  parameters  to  modify  or  supplement  the  altitude 

profiles  of  temperature  and  pressure,  water  vapor, 
ozone,  CH4,  N20,C0 
:  default  was  used  (0) 

-  use  default  profile  for  CO2,  O2,  NO,  SO2,  NO2,  NH3,  HNO3 
:  default  was  used(O) 

-  additional  switch  for  user  defined  atmosphere 
;  set  to  normal  operation  (0) 

-  switch  to  minimize  printing  of  transmittance,  radiance  and  atmospheric 
profile  tables 

;  set  to  minimise  (1) 

-  boundary  temperature  in  K 

;  set  to  sea  temperature.  Note  :  the  sea  temperature  is  only  used  if  the  path 
intercepts  the  earth 

-  surface  albedo  of  the  earth  at  the  location  and  average  frequency  of  the 
calculation  ( 0.0  to  1.0) 

.•  set  to  0 

CARD  2  (line  2)  -  Format  (615,  5F 1 0.3) 

IHAZE,  ISEASN,  IVULCN,  ISCSTL,  ICLD,  IVSA,  VIS,  WSS,  WHH,  RAINRT, 
GNDALT 

IHAZE 

ISEASN 

IVULCN 


-  selects  the  aerosol  model  for  extinction  calculation  and  default 
meteorological  range 

.•  Navy  aerosol  model  used  (3) 

-  selects  the  seasonal  dependence  of  tropospheric  and  stratospheric  profile 
.•  set  to  default  (0)  i.e.  model  selected  in  card  1 

-  selects  the  aerosol  profile  above  the  stratosphere  to  100  km 
;  default  stratospheric  background  used  (0) 


MDEF 

IM 

NORPT 

TBOUND 

SALE 


70 


ICSTL  -  air  mass  character  used  only  with  Navy  aerosol  model  (1  for  open  ocean 

to  10  for  strong  continental  influence) 

;  calculated  values  from  radon  gas  concentration  used  (see  Appendix  F) 
ICLD  -  specifies  cloud  model  and  rain  rates 

;  no  clouds  or  rain  (0)  used 

IVS  A  -  determines  the  use  of  Army  vertical  structure  algorithm  for  aerosols  in 

boundary  layer 
;  not  used  (0) 

VIS  -  specifies  meteorological  range 

;  default  used  (0) 

WSS  -  specifies  current  wind  speed  in  m/s 

;  meteorological  data  used  (see  Appendix  F) 

WHH  -  24  hour  average  wind  speed  in  m/s 

;  meteorological  data  used 
RAINRT  -  rain  rate  in  mm/hr 

;  default  used  (0) 

GND  ALT  -  altitude  of  surface  relative  to  sea  level  in  km 

;  set  to  0 

CARD  3  (line  3)  -  Format  (6F10.3, 15) 

HI,  H2,  ANGLE,  RANGE,  BETA,  RO,  LEN 

HI  -  initial  altitude  in  km  (in  radiance  mode  this  always  defines  the  altitude  of 

the  sensor) 

;  0. 01 7  km  for  pixel-sensor  path  and  pixel  height  for  sea  and  sky  path 
H2  -  final  altitude  in  km 

;  pixel  height  for  pixel-sensor  path  and  0  for  sea  and  sky  path 
ANGLE  -  zenith  angle  measured  from  HI  in  degrees  (between  0  to  1 80),  0  degrees 
is  looking  straight  up 
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;  calculated  for  each  pixel 
RANGE  -  path  length  in  km 

;  not  needed  as  the  first  parameter  is  sufficient  to  specify  the  geometry 
BETA  -  earth  center  angle  subtended  by  HI  and  H2  in  degrees 

;  not  needed  as  the  first  parameter  is  sufficient  to  specify  the  geometry 
RO  -  radius  of  the  earth  at  the  particular  geographical  location  in  km 

;  default  radius  of  earth  i.e.  6371.23  km  (0)  used 
LEN  -  select  normal  or  long  path  through  tangent  height  calculations 

;  normal  operation  (0)  used 

PSI  -  azimuth  of  wind  direction  in  degrees  measured  positive  East  of  North 

from  observer  line  of  sight 
.•  meteorological  data  used 

Seaswitch  -  switch  for  turning  sea  modification  on  or  off 
;  set  to  True  (T) 

CARD  3A1  (line  4)  -  Format  (8F10.3) 

(required  when  lEMSCT  =2  ) 

IPARM,  IPH,  IDAY,  ISOURC 

IP  ARM  -  controls  the  method  of  specifying  solar/lunar  geometry  on  CARD  3A2 

;  observer  lot.  and  long,  used  (1) 

IPH  -  determines  the  type  of  phase  function  to  be  used 

;  MIE  generated  database  for  LOWTRAN  models  (2)  used 
ID  AY  -  day  of  the  year  in  GMT 

;  calculated from  meteorological  data 
ISOURC  -  type  of  extra-terrestrial  source  to  be  considered 
;  sun  (0)  selected 
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CARD  3 A2  (line  5)  -  Format  (8F  10.3) 


PARMl,  PARM2,  PARM3,  PARM4,  TIME,  PSIPO,  ANGLEM,  G 


PARMl 

PARM2 

PARM3 

PARM4 

TIME 

PSIPO 

ANGLEM 

G 


-  observer  latitude  in  degrees  (-90  to  +90;  +  means  North) 

;  32.917  degrees 

-  observer  longitude  in  degrees  (0  to  360) 

:  117.517  degrees 

-  not  used 

-  not  used 

-  GMT  time  in  decimal  hours 

.•  calculated from  meteorological  data 

-  path  azimuth  in  degrees  East  of  North 
.•  calculated from  meteorological  data 

-  phase  angle  of  the  moon 
.•  not  required 

-  asymmetry  factor  for  use  with  HENYEY-GREENSTEIN  phase  function 
.•  not  required 


CARD  4  (line  6)  -  Format  (3F  10.3) 

IVl ,  IV2,  IDV,  IFWHM,  IFILTER 

IV 1  -  initial  frequency  in  wave  number  in  cm"' 

.•  830  cm'^  i.e.  12  used 
IV2  -  final  wave  number  in  cm'^ 

:  1250  cm"’  i.e.  8  /rni  used 
IDV  -  frequency  increment  in  cm'* 

.•  10  cm''  used 

IFWHM  -  frequency  band  at  full  width  half  maximum 
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;  5  cm^  used 

IFILTER  -  specifies  the  filter  to  be  applied 
;  no  filter  (0)  used 


CARD  5  (line  7)  -  Format  (15) 

IRPT 

IRPT  -  determines  what  program  will  do  next 

;  set  to  end  program  (0) 


A  sample  file  is  shown  below. 


F  2  3  2.1 

3  0  0  5 

0.017  0.001 

1  2  100  0 

32.660  117.240 
830  1250 

0 


0  0  0  0 

0  0  0.000 

90.531  0.000 

0.000  0.000 

10  5 


0  0  0  0 

7.000  5.200 

0.000  0.000 

21.181  245.016 

0 


1  287.95  0.00 
0.000  0.000 
0  50.984  T 

0.000  0.000 
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APPENDIX  K.  SEARAD  OUTPUT 


The  format  of  the  SeaRad  output  files  depends  on  the  type  of  calculation  that  is 
done.  Although  all  cases  were  executed  in  thermal  plus  solar  radiance  mode,  their  output 
is  slightly  different  because  of  the  way  the  path  geometry  was  specified.  A  sample  output 
of  each  case  is  below. 


A.  OUTPUT  FOR  AVERAGE  SEA  RADIANCE  PATH 


*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 

DATE:  02/08/1999  TIME:  15:09:32.87 

THERMAL  PLUS  SOLAR  RADIANCE  MODE 

MULTIPLE  SCATTERING  USED 

MARINE  AEROSOL  MODEL  USED 

WIND  SPEED  =  6.00  M/SEC 

WIND  SPEED  =  3.75  M/SEC,  24  HR  AVERAGE 

RELATIVE  HUMIDITY  =  76.11  PERCENT 

AIRMASS  CHARACTER  =  2.0 

VISIBILITY  =  .00  BCM 

SLANT  PATH  TO  SPACE 

HI  =  .011  KM 

HMIN  =  .000  KM 

ANGLE=  118.1001  DEG 

FREQUENCY  RANGE 

IVl  =  830CM-1  (  12.05  MICROMETERS) 

IV2  =  1250  CM-1  (  8.00  MICROMETERS) 

IDV  =  10  CM-1 

IFWHM  =  5  CM-1 

IFILTER=  0 
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SUMMARY  OF  THE  GEOMETRY  CALCULATION 


HI 

H2 

ANGLE  = 
RANGE  = 
BETA  = 
PHI 

HMIN  = 
BENDING  = 
LEN  = 


.011  KM 
.000  KM 
118.100  DEG 
0.023  KM 
.000  DEG 
61.900  DEG 
.000  KM 
.000  DEG 

0 


SEA  AT  287.45  K  REPLACES  BLACK  BODY  BOUNDARY 


UPWIND  =  -5.006  DEG  EAST  OF  LINE  OF  SIGHT 

SUMMARY  OF  OBSERVATION  GEOMETRY 


BETA 

PATH  AZIMUTH 
WIND  AZIMUTH 


00018  DEG 

=  3 14.994  DEG  EAST  OF  NORTH 
=  .  309.988  DEG  EAST  OF  NORTH 


RECEIVER  LATITUDE 
RECEIVER  LONGITUDE 
FOOTPRINT  LATITUDE 
FOOTPRINT  LONGITUDE 
SUBSOLAR  LATITUDE 
SUBSOLAR  LONGITUDE 


=  32.917  NORTH  OF  EQUATOR 
=  117.517  WEST  OF  GREENWICH 
=  32.917NORTH  OF  EQUATOR 
=  117.517  WEST  OF  GREENWICH 
=  7.763  DEG  NORTH  OF  EQUATOR 

=  1 79.866  DEG  WEST  OF  GREENWICH 


VALUES  SEEN  FROM  FOOTPRINT 


10) 


RECEIVER  ZENITH  ANGLE= 
RECEIVER  AZIMUTH 
SOLAR  ZENITH  ANGLE 
SOLAR  AZIMUTH 
SOLAR  SPECULAR  TILT  = 


61.900  DEG 

174.994  DEG  WEST  OF  UP  WIND 
62.650  DEG 

48.828  DEG  WEST  OF  UP  WIND 

40.741  DEG  (  7.03  SIGMA,  PROB  =  1.720E- 


ZERO  RANGE  UNPOLARIZED  VALUES 


SEA  EMISSION  =  29.68807  W  M-2  SR-1  (AV.  EMISS.  .9426) 
SKY  REFLECTION  =  1.05408  W  M-2  SR-1 
SUN  GLINT  =  .00000  W  M-2  SR- 1 
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TOTAL  RADIANCE  =  30.74216  W  M-2  SR-1 
BLACK  BODY  TEMP.=  13.0  C 

FULL  RANGE  UNPOLARIZED  VALUES 

SEA  EMISSION  =  29.14501  W  M-2  SR-1 

SKY  REFLECTION  =  1.03208  W  M-2  SR-1 

SUN  GLINT  =  .00000  W  M-2  SR-1 

PATH  TO  FOOTPRINT  =  .64978  W  M-2  SR-1  (AV.  TRANS.  .9792) 

TOTAL  RADIANCE  =  30.82687  W  M-2  SR-1 
BLACK  BODY  TEMP.  =  13.2  C 


HORIZONTAL 

VERTICAL 

(H-V)/(H+V) 

(WM-2  SR-1) 

(W  M-2  SR-1) 

(%) 

ZERO  RANGE  POLARIZED  VALUES 

SEA  EMISSION 

14.26001 

15.42806 

-3.9 

SKY  REFLECTION 

.85931 

.17277 

66.5 

SUN  GLINT 

.00000 

.00000 

.0 

TOTAL  RADIANCE 

15.13789 

15.60427 

-1.5 

BLACK  BODY  TEMP.  (C) 

-22.0 

-20.7 

FULL  RANGE  POLARIZED  VALUES 

SEA  EMISSION 

14.00002 

15.14499 

-3.9 

SKY  REFLECTION 

3.10762 

2.22232 

16.6 

SUN  GLINT 

.00000 

.00000 

.0 

PATH  TO  FOOTPRINT 

.32489 

.32489 

.0 

TOTAL  RADIANCE 

15.18422 

15.64265 

-1.5 

BLACK  BODY  TEMP.  (C) 

-21.9 

-20.6 

B.  OUTPUT  FOR  PIXEL  TO  SENSOR  PATH 


*****  SEARAD,  A  MODIFICATION  OF  LOWTRAN7  ***** 

DATE:  02/08/1999  TIME:  15:09:10.40 
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THERMAL  PLUS  SOLAR  RADIANCE  MODE 

MULTIPLE  SCATTERING  USED 

MARINE  AEROSOL  MODEL  USED 

WIND  SPEED  =  6.00  M/SEC 

WIND  SPEED  =  3 .75  M/SEC,  24  HR  AVERAGE 

RELATIVE  HUMIDITY  =  76. 1 1  PERCENT 

AIRMASS  CHARACTER  =  2.0 

VISIBILITY  =  .00  KM 

SLANT  PATH,  HI  TO  H2 

HI  =  .017  KM 

H2  =  .011  KM 

ANGLE  =  90.367  DEG 

RANGE  =  .000  KM 

BETA  =  .000  DEG 

LEN  =  0 

FREQUENCY  RANGE 

IVl  =  830  CM-1  (  12.05  MICROMETERS) 

IV2  =  1250  CM-1  (  8.00  MICROMETERS) 

IDV  =  10  CM-1 

IFWHM  =  5  CM-1 

IFILTER  =  0 

SUMMARY  OF  THE  GEOMETRY  CALCULATION 

HI  =  .017  KM 

H2  =  .011  KM 

ANGLE  =  90.363  DEG 

RANGE  =  .948  KM 

BETA  =  .009  DEG 

PHI  =  89.642  DEG 

HMIN  =  .011  KM 

BENDING  =  .000  DEG 

LEN  =  0 

TBOUND  SET  TO  .  1 0  K  FOR  MARINE  SKY 

INTEGRATED  ABSORPTION  =  108.01  CM-1  FROM  830  TO  1250  CM-1 
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AVERAGE  TRANSMITTANCE  =  .7428 


MAXIMUM  RADIANCE  =  3.722E-02  W  M-2  SR-1  (CM-l)-l  AT  830.0  CM-1 

MINIMUM  RADIANCE  =  1.351E-02  W  M-2  SR-1  (CM-l)-l  AT  1090.0  CM-1 

BOUNDARY  TEMPERATURE  =  .10  K 

BOUNDARY  EMISSIVITY  =  1.000 

FILTERED  RADIANCE  =  8.751E+00  W  M-2  SR-1 

BLACKBODY  TEMPERATURE  =  -44.0  C 
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APPENDIX  L.  STORAGE  FORMAT  FOR  SKIN  TEMPERATURE 

DATA 


The  skin  temperature  of  the  ship  was  recorded  through  14  thermistors  attached  to 
selected  locations  around  RA^  POINT  SUR.  The  temperatures  are  read  every  20  seconds 
and  the  original  files  are  identified  by  the  filename  "flirlS  /*Apr"  where  *  is  the  day. 

The  starting  time  and  day  of  the  recording  is  00:00  local  time  (PST)  on  the  first 
day  of  each  month.  At  top  each  file,  the  starting  time  and  day  of  the  first  line  of  record  of 
the  file  is  also  indicated.  For  each  line  of  record,  there  are  17  entries  separated  by  spaces 
across  a  line  of  text.  At  the  end  of  the  record,  a  carriage  return  and  line  feed  character 
starts  the  next  record  on  a  new  line.  The  first  entry  is  the  local  time  in  seconds  from  the 
beginning  of  the  first  day  of  the  month,  and  this  is  followed  by  the  readings  of  the  16 
thermistors  in  ascending  order.  Readings  1  and  9  are  reference  readings  and  should 
always  be  at  25.0.  All  temperatures  are  in  degrees  Celsius. 

A  sample  of  the  data  is  below  : 

april9.1 

Mon,Day,Year=  4  9  1996 
Hour,Min,Sec=  20  8  27 
Volts,istart,istop,nsam,iref  1.49939  1  16  100  1 

763707  25.001  14.458  17.209  13.710  14.943  16.570  13.946 

14.468  25.001  15.562  14.130  31.901  15.310  27.986  15.783 

15.153 
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APPENDIX  N.  ORIGINAL  FILES  AND  R/L  SETTINGS 


Filename 

Original  filename 

Range  setting  (R) 

Level  setting  (L) 

BaseOl.pte 

NSCA0912.ptw 

5 

48 

NSCA0913.ptw 

10 

43 

Base02.pte 

NSCA0915.ptw 

5 

49 

NSCA0916.ptw 

10 

41 

BaseOS.pte 

NSCA0918.ptw 

5 

49 

NSCA0919.ptw 

10 

43 

Base04.pte 

NSCA0927.ptw 

5 

49 

NSCA0928.ptw 

10 

43 

BaseOS.pte 

NSCA0930.ptw 

10 

48 

NSCA0931.ptw 

10 

43 

Base06.pte 

NSCA0933.ptw 

5 

49 

NSCA0934.ptw 

10 

43 

BaseOT.pte 

NSCA0936.ptw 

5 

49 

NSCA0937.ptw 

10 

43 

BaseOS.pte 

NSCA0938.ptw 

5 

47 

NSCA0939.ptw 

10 

43 

Base09.pte 

NSCA0940.ptw 

5 

47 

NSCA0941.ptw 

10 

44 

BaselO.pte 

NSCA0942.ptw 

5 

47 

NSCA0943.ptw 

10 

43 

Basel  l.pte 

NSCA0944.ptw 

5 

47 

NSCA0945.ptw 

10 

44 

Basel2.pte 

NSCA0946.ptw 

5 

48 

NSCA0947.ptw 

10 

44 

BaselS.pte 

NSCA0948.ptw 

5 

47 

NSCA0949.ptw 

10 

43 
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Filename 

Original  filename 

Range  setting  (R) 

Level  setting  (L) 

Basel4.pte 

NSCA0950.ptw 

5 

47 

NSCA0951.ptw 

10 

44 

BaselS.pte 

NSCA0952.ptw 

5 

46 

NSCA0953.ptw 

10 

43 

Basel6.pte 

NSCA0954.ptw 

5 

46 

NSCA0955.ptw 

10 

44 

Basel  7.pte 

NSCA0956ptw 

5 

47 

NSCA0957.ptw 

10 

44 

Base27.pte 

NSCA0978.ptw 

5 

49 

NSCA0979ptw 

10 

43 

Base28.pte 

NSCA0980.ptw 

10 

51 

NSCA0981.ptw 

10 

44 

Base29.pte 

NSCA0982.ptw 

5 

51 

NSCA0983.ptw 

10 

45 

BaseSO.pte 

NSCA0988.ptw 

5 

51 

NSCA0989.ptw 

10 

44 

Base3 1  .pte 

NSCA0990.ptw 

5 

47 

NSCA0991.ptw 

10 

43 

Base32.pte 

NSCA0993.ptw 

5 

49 

NSCA0994.ptw 

10 

43 

Base33.pte 

NSCA0996.ptw 

5 

49 

NSCA0997.ptw 

10 

43 

Base34.pte 

NSCA0999.ptw 

5 

49 

NSCAlOOO.ptw 

10 

43 

Base35.pte 

NSCA902.ptw 

5 

49 

NSCA903.ptw 

10 

43 

Base36.pte 

NSCA905.ptw 

5 

50 

NSCA906.ptw 

10 

43 
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Filename 

Original  filename 

Range  setting  (R) 

Level  setting  (L) 

Base37.pte 

NSCA911.ptw 

5 

50 

NSCA912.ptw 

10 

44 

Base38.pte 

NSCA914.ptw 

5 

50 

NSCA915.ptw 

10 

43 

Base39.pte 

NSCA917.ptw 

10 

51 

NSCA918.ptw 

10 

43 

Base40.pte 

NSCA920.ptw 

10 

51 

NSCA921.ptw 

10 

44 

Base41.pte 

NSCA923.ptw 

5 

50 

NSCA924.ptw 

10 

42 

Base43.pte 

NSCA929.ptw 

5 

49 

NSCA930.ptw 

20 

44 

Base44.pte 

NSCA932.ptw 

5 

49 

.  NSCA933.ptw 

20 

43 

Base45.pte 

NSCA935.ptw 

5 

49 

NSCA936.ptw 

20 

47 

Base49.pte 

NSCA947.ptw 

5 

49 

NSCA948.ptw 

10 

43 

BaseSO.pte 

NSCA958.ptw 

5 

51 

NSCA959.ptw 

10 

44 
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APPENDIX  O.  AN  EXTRACT  ON  PERFORMING  COMPLEX 
THERMAL  MEASUREMENTS  FROM  THE  AGA  MANUAL  [Ref.9] 


Section  10  of  the  AGA  operating  manual  describes  the  thermal  measurement 
techniques  for  making  measurements  with  the  AGA  Thermovision  infrared  scanner.  The 
extract  in  the  next  three  pages  describes  the  theory  and  how  complex  thermal 
measurements  are  made  with  the  AGA  scanner. 
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area)  is  jUSt  brightened  up  by  :he  isothernrj.  For 
ease  of  operation  this  is  often  earned  out  with  the 
PICTURE  MODE  switch  set  to  INVERTED  or  BLACK. 

(c)  Adjust  the  second  isotherm  (LEVEL  2)  to  indicate 
the  colder  of  the  two  areas  of  interest.  If  this  cannot 
be  achieved,  increase  the  THERMAL  RANGE  one 
step  at  a  rime  and  repeat  (a)  and  (b)  and  (c)  again 
until  both  colder  and  warmer  areas  are  covered  by 
the  same  selected  range. 

NOTE;  For  overall  accuracy  maintain  the  THERMAL 
RANGE  setting  at  its  lowest  level  consistent  with 
displaying  both  the  reference  and  object  area 
isotherms. 

Without  readjusting  THERMAL  RANGE  or  LEVEL 
proceed  as  follows: 

\d)  Carefully  adjust  the  two  isotherm  controls  to 
indicate  the  object  area  to  be  measured  and  the 
reference  area  respectively.  This  is  shown  by  the 
areas  being  displayed  as  saturated  white. 

(e)  Read  the  isotherm  scale  values  at  the  top  edge  of 
both  isotherm  markers,  calculate  the  difference 
and  multiply  by  the  THERMAL  RANGE  setting  (see 
figure  2). 

This  results  in  a  thermal  difference  value  that  can 
be  applied  directly  to  the  calibration  curve  since  we 
have  postulated  e  and  ^^-1.  Proceed  as 

follows: 

if)  Mark  the  reference  temperature  on  the  calibration 
cnart.  Read  via  the  curve  the  corresponding 
thermal  value  ij..  Add  the  thermal  difference  value 
—Iq,.  obtained  in  (e)  to  l^.  This  gives  the  blackbody 
thermal  value  Iq  of  the  object,  which  is  then 
converted  to  temperature  t^  by  means  of  the  curve. 

NOTE:  If  the  reference  should  have  higher  temperature 
than  the  object,  the  thermal  difference 
becomes  negative  and  should  consequently  oe 
subtracted  from  the  reference  thermal  value  The 
lower  half  of  the  isotherm  scale  has  negative  readings, 
it  is  therefore  important  to  note  the  correct  sign  of  the 
isotherm  scale  readings  when  calculating  -  v. 

COMPLEX  THERMAL  MEASUREMENTS 

The  previous  instructions  have  assumed  the  ideal 
situation,  that  the  object  is  a  perfect  blackbocy  and  mat 
no  external  factors  influence  the  measurements. 

In  complex  thermal  measurement  situations,  however, 
where  several  factors  influence  the  measurements,  me 
true  object  temperature  has  to  be  derived  by 
calculation. 

An  expression  is  set  up  to  include  the  radiation  that 
reaches  the  scanner  from  the  ooject.  From  this 
expression  the  unknown  parameter  is  extracted  ;e.g. 
object  temperature  or  object  emissivityi.  Before  mis 
derivation  is  carried  out  the  following  variables  should 
bev^onsidered. 

Variables  affecting  measurement 

The  most  important  external  influencing  factors  are: 
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Emissivfty  -  Actual  objects  are  seldom  "biack"  The 
emissivity  factor  must  therefore  often  be  taken  into 
account  m  infrared  temperature  measurements  as 
detailed  m  the  measurement  formulae.  Individual 
object  emissivity  can  be  measured  or  can  be  found  m 
tables,  (e.g.  Table  1  page  n.3). 

Ordinarily,  object  materials  and  surface  treatments 
exhibit  emissivities  ranging  from  approximately  0.1  to 
0-95.  A  highly  polished  (mirror)  surface  falls  below  0. 
while  an  oxidised  or  painted  surface  has  greatly 
increased  emissivity.  Oil-based  paint,  regardless  of 
colour  in  the  visible  spectrum,  has  an  emissivity  over 
0.9  in  the  infrared.  Human  skin  exhibits  an  emissivity 
close  to  1  which  is  used  of  in  medical  thermography  for 
direct  body-temperature  readings. 

Surroundings  -  The  reflectivity  factor  of  an  epaoue. 
diffuse  surface  is  ■=  1  -  e  . 

A  low  emissivity  factor  means  therefore  not  only  mat 
the  emission  from  the  object  is  lower  than  from  a 
blackbody  of  the  same  temperature.  It  also  means  that 
undesired  radiation  from  the  surroundings  is  reflected 
into  the  scanner.  The  measurement  formulae  ’S 
corrected-  for  reflected  radiation  by  means  of  a  term 
including  I3,  which  is  related  to  ambient  temperature,  it 
IS  therefore  important  that  the  object  surroundings 
have  a  homogenous  (ambient)  temperature  and  do  not 
include  hot  areas  so  positioned  that  their  radiation  can 
be  reflected  by  the  object.  Sometimes  an  efficient 
ambient  temperature  has  to  be  estimated  :o  take  into 
account  radiation  sources  that  cannot  be  removed  or 
shielded. 

Object  opacity  -  Objects  may  be  more  or  less 
transparent  to  infrared  radiation  e.g.  glass  anc  plasncs. 
The  .''adiation  received  from  such  an  object  would 
include  irrelevant  transmitted  radiation  from  the 
background.  Besides  the  efficient  emissivity  would  be 
low  and  hard  to  estimate.  However,  the  spectral 
cnaracteristics  of  the  scanner  may  be  adapted  to  tne 
object  by  means  of  spectral  filtering  such  that  only 
spectrum  regions  where  the  emissivity  is  hign  and  the 
transmission  is  zero  are  utilised.  Where  this  is  possible 
to  do,  thermal  measurements  can  be  carried  out. 
However,  the  resultant  emissivity  and  reflectivity  of  tne 
object  must  still  be  taken  into  account  if  the  spectral 
filtering  does  not  result  in  an  effective  blackbody. 

Object  size  -  The  detector  suotends  a  certain  solid 
angle  in  the  object  space.  !f  the  object  does  not  cover 
this  angle  the  detector  will  .'•eceive  radiation  from  the 
object  background.  This  will  tend  to  lower  the 
temperature  difference  indicated  berween  the  cDiect 
and  background.  It  is  therefore  recommended  that  the 
size  of  the  image  area,  where  temperature  s  to  be 
measured,  snould  not  be  less  than  2.5  x  2.5  mm  as  seen 
on  the  Black/White  monitor  chassis  display  screen. 

The  Atmosphere  -  Certain  constituents  of  the  atmos¬ 
phere  absorb  infrared  radiation  in  the  spectral  bands 
being  used.  The  most  important  gases  are  water 
vapour  (H,0}  and  carbon  dioxide  (CO>i.  This  absorption 
will  attenuate  the  infrared  radiation,  from  the  oDject  to 
the  instrument.  In  some  cases  -e.g.  hot  comeustion 


in 
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flames)  the  self-emission  of  atmospheric  gases  is  verv 
high  and  has  to  be  taken  into  consideration. 

The  atmospheric  anenuation  can  be  calculated  and 
compensated  with  great  accuracy  using  either 
computer  programs  or  an  approximative  formula.  iSee 
Section  10  page  10.6.) 


Measurement  formulae 

in  general  a  number  of  different  parameters  are 
involved  in  temperature  determinations  using 
radiometric  methods.  Here  are  listed  all  parameters 
used  in  the  following  description.  They  are  divided  into 
three  groups  of  which  the  last  is  used  only  for  formula 
derivation. 


Radiation  source  and  atmosphere  related  parameters 


^0 


^atm 

e 


'  0 


1 


I. 

^3tm 


object  temperature  t  Ci 

temperature  of  object  surroundings 
(ambient)  ('O 

reference  temperature  i  C) 
atmosphere  temperature  ("O 

emissivity  factor,  0^^  e  'S 1 
(e  =1  for  a  "blackbody  "} 

object  emissivity 
emissivity  of  reference  surface 

emissivity  of  object  surroundings 
6  3  is  considered  ^1 

atmospheric  correction  factor  for  the 
atmosDhere  between  scanner  and  object, 
i.e.  atmospheric  transmission  corrected 
for  systems  spectral  resoonse 

as  above  for  the  atmosphere  between 
scanner  and  reference 

distance  scanner  to  object  (ml 
as  above  to  reference  (m) 

calibrated  thermal  value.  Refers  to  the 
instruments  calibration  function 
^graphical  curve  or  calculator  program; 
which  describes  the  relation  between 
object  blackbody  temperature  and  instru¬ 
ment  numerical  output  on  short  distance. 
I=f(t)  (lU) 

calibrated  thermal  value  for  t=tQ  (lU) 
calibrated  thermal  value  fort^tg  ’lU) 
calibrated  thermal  value  for  t=tj.  (lU) 
calibrated  thermal  value  fort  =  t3^Pn  '!U) 


Instrument  reading  parameters 

-  measured  object  thermal  /aiue.  instru¬ 
ment  reading  of  radiation  received  from 
the  object  surface  via  the  atmosphere  ilU) 


ij.'  -  as  above  for  the  reference  (lU) 

The  instrument  reading  consists  of  two  terms:  i'  =L-i, 
where 

L  -  THERMAL  LEVEL  setting  of  the  instru¬ 

ment  (lU) 

Lq  -  as  above  at  object  measurement  dU) 

Lj.  =  as  abov^feference  measurement  (lU) 

i  =  measured  relative  thermal  value  = 

instrument  scale  or  display  reading  x 
THERMAL  RANGE  setting  (lU) 

=  as  above  referring  to  object  measurement 
(lUl 

ip  =  as  above  referring  to  reference  measure¬ 

ment  (lUl 

=  as  above  referring  ro  the  difference 
between  object  and  reference  measure¬ 
ment  lIU) 

— 'or“’o“'r 


Radiation  notations 


received  and  detected  radiation  from  a 
blackbody  of  temperature  t,  at  snort 
distance  i  Photons/ sec) 


$3  =  as  above  when  t=t3  ?Photons/sec; 

^3tm  ”  as  apove  when  t=t3jpp  (Photons/sec) 

=  radiation  received  via  the  atmospnere 
from  the  object  surface  i  Photo  ns/ sec  i 


Rg.  3  Radiation  conditions  in  the  general 
measurement  situation 

.^igure  3  illustrates  schematically  the  radiation 
conditions  in  the  general  measurement  situation.  .Note 
that  the  reflection  factor  of  the  object  surface  is  1  t  q 
for  an  opaque  object,  ana  that  the  emission  factor  of  the 
atmosphere  is  (1  -  ty.  With  the  support  of  this  figure 
the  received  radiation  can  be  written: 

Sq  —  Td  ■  ^  ^ .0^  ^3  ^  ^atm 

=oDjecr  radiation-hreffected  radiation - 
atmosphere  radiation 


10.4 
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This  radiation  relationship  can  be  converted  into  a 
thermal  value  relationship  by  utilisinq  the  fact  that  AGA 
Thermovision'^  780  has  a  linear  ,  ...ton  counting 
detector. 

We  can  therefore  write  the  thermal  value 


in  a  practical  measurement  situation,  we  assume  that 
object  and  reference  are  placed  close  together.  Hence 


r,  =  r  and 


Ur-U-  With  these  assumotions  we 


ao”'ar“'a' 


where  C  is  an  empirical  instrument  factor 
Replacing  Sq  ,Sq,  Sg  and  S3fpp|  by 
respectively  results  in 


This  is  the  general  measurement  formula  describing 
how  the  measured  thermal  value  Iq'  relates  to  the  three 
radiation  sources:  object,  surroundings  and 
atmosphere,  together  with  their  associated  parameters 
emissivity  and  transmission.  The  object  is  considered 
opaque  (i.e.  no  transmitted  radiation  through  the 
object)  and  all  the  relevant  surrounding  surfaces  are 
assumed  to  have  the  same  (ambient)  temperature  t3 
(from  that  assumption  follows  automatically  e3=1) 

The  radiation  terms  are  expressed  as  blackbooy 
thermal  values  determined  by  the  calibration  function 
at  respective  source  temperature. 

This  formula  is  a  suitable  starting  point  for  derivation  of 
special  measurement  formulae  adapted  ro  special 
measurement  situations. 

For  direct  measurements  {without  reference)  the 
formula  can  be  rewritten,  using  the  relation  iQ  =  LQ^iQ. 
Solve  for  IqI 


>  -  *-0**0  .1  ,,,  1,1 

'*0  ^  ^ 

In  case  the  object  temperature  is  known  but  the  ooject 
emissivity  unknown,  the  latter  can  oe  calculated.  Solve 


"o  ■  ‘o  'o  ^ 


r  (I  -  !  ) 

0  0  3 

For  relative  measurements  with  a  known  temperature 
reference,  corresponding  formulae  are  derived  as 
follows: 

!f  =  -r,  C-V'ar  - 

Where  IgQ  and  ‘3^  refer  ro  the  surroundings  of  oDject 
and  reference  respectively. 


*r  ^o^’o  'r  ■“  ‘or 

As  previously  described  under  relative  measurements, 
'the  THERMAL  LEVEL  should  be  keot  constant  wnen 
using  this  method.  Thus  L^^  -' L^  and  from  that  follows: 


■'r  =  ^'or 


.^1,  Ml  --il  I3 

tV  e 


Solve  for  i  qI 


Cr-I.) 


The  above  de rivaled  formulae  can  be  a  great  deal 
simplified  under  special  conditions.  For  instance  if  all  e 
and  r  are  close  to  unity  we  get 

^o  ^0  ^0  ^  “^or 

for  direct  and  relative  measurement  respectively,  i.e.  as 
assumed  in  the  first  part  of  this  section. 

The  formula  for  relative  measurement  will  be  reduced 
and  thereby  the  overall  accuracy  Improved,  if 
and  r^. 

This  requirement  is  met  if  for  example  the  reference 
surface  is  a  part  of  the  object  surface  itself  where  the 
temperature  is  known.  'The  formula  is  then  reduced  to 


ATMOSPHERIC  INaUENCS  ON 
INFRARED  TEMPERATURE  MEASUREMENTS 

The  atmosphere  is  not  fully  transparent  ro  infrared 
radiation.  !t  attenuates  radiation  passing  through  it  and 
also  emits  radiation,  rne  amount  dependent  on  the 
temperature. 

The  measurement  formulas  previously  described 
indicated  now  to  correct  for  the  atmospheric  influence 
provided  rhai  the  correction  factor  r  is  known,  r 
depends  on  a  number  of  parameters.  The  most 
important  parameters  are  listed  below  together  witn 
the  actual  values  for  what  -s  defined  as  "standard 
atmosphere".  iSee  Table  1.) 

Table  1  Parameters  influencing  the 
atmospheric  correction  factor  r 


Parameter  Standard  atmosphere  j 

Scanner  type  (SW,  38AR.  LW) 

Distance  scanner  to  object 
j  Air  Temperature  15C 

I  Air  pressure  larm  =  10l3mD 

Relative  humidity  50% 

Carbon  dioxide  (CO,)  300  ppm 

Carbon  monoxide  ((30  t  ppm 

Oinitrogen  oxide  (NjO)  0.3  ppm 

Ozone  {0.j)  0.355  ppm 

Methane  (CHJ  2  ppm 

Aerosol  visibility  iparticies)  10  km 
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